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I .  INTRODUCTION 


Nitrocellulose-based  propellants  are  widely  used  throughout  the  military 
and  civilian  areas  in  guns  and  rockets.  There  are  almost  as  many  different 
procedures  for  analysis  as  there  are  propellants.  Even  in  the  area  of  liquid 
chromatography,  frequently,  more  than  one  procedure  is  used  in  the  analysis  of 
the  propellant.  Most  techniques  have  been  designed  to  answer  a  specific 
question  about  the  purity  or  amount  of  two  or  three  ingredients.  However, 
most  propellants  have  at  least  three  or  four  original  Ingredients  and  fre¬ 
quently  one  or  two  other  reaction  products  that  can  be  easily  extracted. 

Since  these  ingredients  can  be  extracted  into  a  solvent  it  would  seem 
logical  to  try  to  analyse  for  all  of  them  at  one  time.  Early  efforts  centered 
on  developing  separations  using  silica  columns  such  as  Corasil  II  and  eluting 
with  non  polar  solvents  such  as  hexane,  methylene  chloride,  and  dioxane  [1-14]. 
A  major  difficulty  to  this  in  the  past  has  been  that  silica  columns  are  dif¬ 
ficult  to  work  with.  Precise  control  of  water,  alcohols,  and  other  polar 
ingredients  is  necessary  to  obtain  reproducible  retention  times  and  amounts. 

II.  BACKGROUND 

Chromatographic  efforts  shifted  to  reverse  phase  separations  with  the 
advent  of  permanently  bonded  phases  and  improved  equipment.  Reverse  phase 
columns  do  not  have  the  problems  that  silica  columns  do.  Materials  are 
separated  on  differences  in  their  nonpolar  segments  rather  than  in  their  polar 
segments  a3  in  silica  columns.  Currently  over  80%  of  all  HPLC  separations 
are  performed  on  reverse  phase  columns  [15],  but  investigators,  have  expressed 
difficulty  in  obtaining  good  separations  of  complex  mixtures  with  the  binary 
solvent  mixtures  commonly  used  [16-18].  Nonetheless,  analysis  times  have  been 
shortened  from  10  to  12  hours  using  the  older  MIL-STD-286B  methods  to  less 
than  one  hour  using  HPLC  [19-20]. 

Recently,  significant  interest  has  beer  observed  in  using  more  than  2 
solvent  systems  with  C-18  columns.  Manuf ac . jrers  have  started  producing  HPLC 
equipment  that  can  pump  3  or  4  solvents  simultaneously.  The  reason  for  a 
multi-solvent  system  is  that  separations  are  accomplished  by  more  than  dif¬ 
ferences  in  polarity.  Snyder  [21]  has  proposed  that  solvents  should  be 
classified  as  to  three  properties:  proton  acceptor,  proton  donor,  and  dipole 
interaction.  It  is  believed  that  these  properties  make  the  major  contribu¬ 
tions  to  separations  in  HPLC.  Kirkland  [22]  has  developed  a  procedure  con¬ 
sisting  of  7  experiments  which  can  be  used  to  predict  the  best  solvent 
composition1  for  separating  ingredients  in  a  mixture.  A  program  was  undertaken 
to  develop  a  procedure  that  will  separate  the  maximum  number  of  commonly  used 
propellant  Ingredients  possible. 

Previously  [23]  a  procedure  was  described  where  with  a  series  of  7  experi¬ 
ments  thj  separation  of  up  to  10  propellant  ingredients  would  be  optimized. 

The  procedure  is  not  one  that  should  be  used  for  every  unique  analytical 
problem  encountered  but  shows  great  promise  for  the  occasionally  difficult 
cases  and  for  repetitive  analysis.  In  production  facilities  this  latter  case 
will  be  of  significant  aid.  Frequently,  several  samples  of  the  same  material 
require  rapid  analysis  each  day.  It  is  not  unusual  for  more  than  one  for¬ 
mulation  to  be  in  production  at  the  same  facility.  Isocratic  separations  are 
desired  for  their  reproducibility  and  turn  around  time,  also  no  column 


reequilibration  is  required.  If  one  set  of  conditions  could  be  used  for  all 
the  formulations  manufactured  at  the  facility,  the  samples  could  be  analysed 
in  a  random  or  as  received  order  rather  than  in  groups. 

At  the  time  of  writing  the  first  paper  on  this  subject  several  limita¬ 
tions  still  existed.  The  computer  program  employed  could  only  work  with  10 
ingredients  and  was  limited  to  a  64K  mini  computer  operating  in  Basic.  This 
limitation  was  imposed  by  computer  memory  restrictions.  A  disc  based  computer 
was  recently  made  available  to  the  RD&E  Center  laboratory.  With  the  computers 
increased  capabilities  for  large  matrix  operations  it  was  possible  to  modify 
the  program  to  consider  many  more  materials  simultaneously.  A  limit  of  35 
ingredients  was  arbitrarily  selected  as  the  maximum  number  for  calculations. 

A  program  this  size  performs  a  complete  evaluation  to  595  pairs  of 
ingredients.  As  a  practical  limit  20  propellant  ingredients  were  selected  for 
further  work. 

A  new  disc  based  program  has  now  been  written.  It  includes  several 
improvements  over  the  original  program  such  as  calculating  with  the  actual 
peak  widths  rather  than  an  average,  and  plotting  all  resolution  lines  on  one 
graph  on  a  x-y  plotter.  This  program  has  been  written  in  Basic.  A  Fortran  77 
version  is  also  available.  See  Appendix  B  for  Basic  programs  "XLDB6  and 
XLDB6B",  see  Appendix  C  for  Fortran  77  program  "LC  Main". 

III.  EXPERIMENTAL 

The  chromatographic  system  consisted  of  a  Waters  Associates  720  System 
Controller,  710B  WISP  sample  injector*  three  M6000  pumps,  a  Perkin-Elmer  LC  85 
variable  wavelength  ultra  violet  detector  and  either  a  Supleco  250  mm  x  4.6  mm 
ID  Supelcosil  5  micron  C-18  column  or  a  Dupont  5  micron  ZORBAX  ODS  4.6  mm  x 
250  mm  column.  Solvents  were  from  Burdick  and  Jackson  and  were  filtered  prior 
to  mixing.  Solvents  in  this  study  were  combinations  of  the  identified  solvent 
and  distilled  deionized  water.  MeOH  is  60%  methanol-40%  water,  ACN  is  50% 
acetonitriie-50%  water,  and  THF  is  42%  tetrahydrofuran-50%  water. 

The  analog  UV  signal  was  digitized  and  sent  to  a  Hewlett-Packard  3357 
Laboratory  Automation  System  for  data  reduction.  Following  normal 
integration,  a  computer  program  was  automatically  run  that  calculated  the  peak 
widths  and  retention  times  >*?-’uired  for  input  to  the  solvent  selection  program 
(see  Basic  program  "Peak  WT"  in  Appendix  A). 

Samples  were  prepared  by  cutting  a  propellant  into  small  pieces  and 
accurately  weighing  50-100  ug  into  a  15  ml  test  tube.  An  accurately  weighed 
internal  standard  was  added  to  the  test  tube  and  then  10  ml  of  ,i  50% 
acetonitrile/water  solution  was  added.  An  ultrasonic  micro  tip  was  placed  in 
the  test  tube  and  the  sample  was  extracted  for  15  minutes.  The  solution  was 
filtered  and  placed  in  a  WI0P  cample  vile  for  analysis. 

n  50%  acetonitrile  solution  was  selected  because  noncrcss  linked  nitro¬ 
cellulose  samples  will  dissolve  completely  or  form  a  gum  in  pure  acetonitrile. 
If  this  solution  is  filtered  or  injected  onto  a  column  the  nitrocellulose 
would  plug  the  pores. 
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When  HMX  is  suspected  in  the  sample  and  quantitative  analysis  is 
required,  a  smaller  sample  or  greater  amount  of  solvent  is  needed.  HMX  has 
only  a  limited  solubility  in  even  pure  acetonitrile.  RDX  is  much  more  soluble 
than  HMX  and  should  not  present  a  problem.  You  must  allow  for  the  presence  of 
HMX  and  RDX  even  when  no  qualitative  data  is  required,  when  selecting  the 
optimum  chromatographic  solvent. 

IV.  SOLVENT  SELECTION 

A  column  is  first  characterized  by  determining  the  retention  times  and 
peak  widths  of  all  known  propellant  ingredients  in  7  predetermined  solvent 
compositions.  This  data  is  fed  into  a  computer  file  for  later  reference. 

This  only  needs  to  be  done  once  for  each  column  (see  Table  1).  This  procedure 
has  been  described  earlier  [22]. 

Most  samples  are  submitted  for  analysis  in  RD5.E  Center  laboratory  as  part 
of  an  aging  surveillance  program  or  for  verification  of  composition.  In  these 
cases  the  computer  is  asked  to  analyse  the  data  file  for  only  those  ingre¬ 
dients  known  or  suspected  to  be  present  in  the  propellant  extract.  A  typical 
example  of  this  procedure  is  the  N-5  propellant. 

The  frest  propellant  is  known  to  contain  Nitroglycerine, 

Diethylphthalate,  and  2-Nitrodiphenylamine .  The  binder,  lead  salts,  and  other 
insoluble  Ingredients  are  not  considered.  As  the  propellant  ages,  degradation 
products  of  2-NDPA  are  expected  to  be  produced.  Therefore,  the  computer  was 
asked  to  predict  the  best  conditions  for  separation  of  2,2-Dinitro  DPA, 
2j,4-Dinitro  DPA,  and  2,4,4-Trinitro  DPA,  an  internal  standard 
Dipropyiphthalate,  as  well  as  the  three  original  ingredients.  It  predicted 
the  best  solvent  composition  for  the  separation  would  be  82%  MeOH/18%  ACN  with 
a  minimum  resolution  of  2.855  for  any  two  peaks  (see  Table  2).  The  triangular 
graph  (Figure  1)  indicated  good  resolution  could  be  obtained  with  MEOH/ACN 
combinations  between  86/14  and  78/22.  Since  methanol  has  a  higher  viscosity 
than  acetonitrile  the  78/22  mixture  was  selected.  This  permitted  a  higher 
flow  rate  than  possible  with  the  82/18  mixture. 

As  can  be  seen  in  the  chromatograph  of  the  standard  sample  (Figure  2), 
all  the  ingredients  were  veil  resolved  with  the  78/22  mixture  at  2.0  ml/min. 
Once  the  system  was  calibrated  the  propellant  samples  were  analysed. 

V.  DISCUSSION 

One  aspect  cf  the  procedure  that  has  been  receiving  attention  is  its  abi¬ 
lity  to  predict  optimized  conditions  for  any  combination  of  ingredients  in  a 
data  base.  Once  the  retention  times  of  a  large  number  of  ingredients  have 
been  determined  for  the  7  solvent  compositions,  the  optimum  conditions  for  any 
combination  of  those  Ingredients  can  quickly  be  predicted  without  further 
experimentation.  A  drawback  to  this  is  the  data  base  for  a  Waters  10  micro¬ 
meter  C-18  radial  pack  is  not  necessarily  compatible  to  another  column  such  as 
a  Waters  10  micrometer  steel  column.  To  explore  how  significant  column  dif¬ 
ferences  effected  the  predicted  optimum  conditions  a  data  base  was  established 
for  several  different  columns.  Tables  3  thru  7  list  the  retention  times  for 
20  propellant  ingredients  in  the  7  solvent  compositions.  As  reported  earlier, 
MeOH  signifies  60%  methanol  in  water,  ACN  signifies  50%  acetonitrile  in  water, 
and  THF  signifies  42%  tetrahydrofuran  in  water.  Table  8  lists  the  important 
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specifications  for  the  5  columns  tested.  Three  Waters  10  micrometer  018 
Radial  Pac  columns  were  tested.  The  computer  predicted  almost  exactly  the 
same  optimum  conditions  for  all  three.  The  obvious  implication  is  that  for 
these  bulk  loaded  columns  the  column  to  column  reproducibility  is  very  good. 

It  should  be  noted  however  that  this  data  was  obtained  using  a  RCM-100  modual. 
The  Z-modual  uses  slightly  different  columns  and  produce  different  results 
even  though  they  have  the  same  name. 

Following  the  procedure  described  earlier,  the  computer  predicted  optimum 
solvent  compositions  were  determined  for  each  of  the  5  different  columns.  As 
seen  in  Figure  3  most  of  the  columns  have  a  optimum  composition  in  the  same 
general  area. 

Not  unexpectedly,  a  significant  improvement  in  resolution  was  obtained  in 
changing  from  a  10  micrometer  packing  to  a  5  micrometer  packing.  Figure  A  is 
a  chromatographic  run  using  a  5  micrometer  Waters  C-18  Radial  Pac  column  at 
the  optimum  solvent  composition.  Figure  5  is  a  run  of  the  same  solvent  but  on 
a  10  micrometer  C-18  Radial  Pac  column. 

Eighteen  ingredients  can  be  separated  well  enough  to  quantify  at  the  pre¬ 
dicted  composition  for  the  5  micrometer  Waters  column  (Figure  A).  Two  pairs 
of  these  peaks  appear  as  fused  peaks  while  the  remaining  ingredients  are  well 
resolved.  Only  ethylcentrallte  and  diphenylamine  are  unresolved  at  this 
composition.  In  Figure  5,  the  10  micrometer  Waters  column  shows  almost  no 
completely  resolved  peaks  for  the  same  conditions. 

At  the  predicted  optimized  conditions  for  the  Perkin-Elmer  HS-5  C-18 
column  we  see  that  several  of  the  later  peaks  have  shifted  (Figure  6). 
Nitroglycerine  has  fused  with  2, A  DNT  and  ethylcentrallte  with  TNT. 
Diphenylamine  has  eluted  later  and  is  almost  resolved  completely  from  BTTN. 
Four  pairs  are  completely  unresolved:  .  resorcinol  and  triacetin;  RDX  and 
methylcentrallte;  2, A  DNT  and  nitroglycerine;  TNT  and  ethylcentrallte. 
Generally  the  ingredients  appear  to  be  less  resolved  than  in  the  Waters  5 
micrometer  radial  pac  system. 

A  suplecosil  LC-18  column  (Figure  7)  was  optimized  for  separating  the 
mixture.  Again  the  predicted  solvent  composition  was  similar  to  the  previous 
examples.  Comparison  with  the  Perkin-Elmer  column  shows  that  ethylcentrallte 
and  TMETN  have  shifted  relative  to  the  other  ingedients.  Diphenylamine  is 
completely  resolved  as  is  TNT  but  now  the  ethylcentrallte  is  co-eluting  with 
BTTN.  The  resolution  of  fused  peaks  is  better  than  in  the  Perkin-Elmer  column 
but  there  are  still  A  pairs  of  co-eluting  peaks. 

The  Dupont  Zorbax  ODS  column  had  a  predicted  optimum  significantly  dif¬ 
ferent  from  the  others.  Figure  8  indicates  complete  resoltuion  of  all  the 
ingredients.  Dupont  literature  discusses  the  unique  properties  of  their 
Zorbax  columns..  They  ara  manufactured  with  monofunctional  silane  units  to 
produce  a  "monolayer" .  Whether  or  not  the  characteristics  of  this  column  are 
unique  or  common  to  all  Zorbax  ODS  columns  will  have  to  be  determined.  The 
long  elution  times  of  the  last  A  ingredients  is  undesirable  but  if  the  analyst 
wishes  to  use  an  isocratic  system  this  situation  is  unavoidable.  A  change  in 
the  solvent  composition  after  17  minutes  could  cause  these  ingredients  to 
elute  fa  .ter  with  minimal  loss  of  resolution;  This  would  introduce  an 
equilibration  delay  between  runs  however,  so  a  trade  off  would  have  to  be 


made.  Using  the  solvent  composition  chosen  for  the  Supelco  column,  this 
Dupont  column  produced  a  chromatogram  similar  to  Figure  7  except  that 
nitroglycerine  and  ethylcentralite  appear  as  fused  peaks  to  2,4  DNT  and  BTTN 
respectively. 

Figure  9  represents  a  chromatogram  of  the  Waters  10  micrometer  Radial  Pac 
column  at  Its  best  composition  near  the  selected  Dupont  composition.  There  is 
some  shifting  of  peaks  but  the  chromatograms  appear  similar.  The  most  impor¬ 
tant  difference  is  the  significant  loss  of  resolution.  There  are  4  fused  peak 
pairs  and  4  sets  of  unresolved  peaks.  While  on  first  glance  Figure  9  appears 
to  be  better  than  Figure  4,  close  examination  indicates  more  partially 
resolved  Ingredients  than  in  Figure  9.  So  the  original  predicted  optimum  com¬ 
position  for  the  Waters  column  is  indeed  better  than  the  alternate. 

VI.  ANALYSIS  OF  UNKNOWN 

Occasionally,  a  propellant  of  unknown  composition  is  submitted  for 
analysis.  A  variety  of  analytical  tools  are  employed  for  this  task. 

For  LC  analysis,  the  solvent  capable  of  separating,  at  least  partially, 
all  known  ingredients  is  initially  selected.  Retention  times  of  the  unknown 
peaks  are  compared  with  the  predicted  retention  times.  The  retention  times  of 
the  authentic  ingredients  are  then  obtained  for  comparison.  Oroe  a  list  of 
probable  Ingredients  is  obtained,  the  computer  selects  a  second  solvent  to 
separate  these  specific  materials.  Standards  are  then  prepared  and  the  analy¬ 
sis  proceeds  normally,  • 

One  foreign  propellant  (I?  1)  analysed  in  this  laboratory  was  found  to  be  a 
rather  simple  double-base  composition.  The  extractables  were  identified  as 
nitroglycerine,  NDPA,  and  raethylcentralite.  Excellent  resolution  was  easily 
obtained  with  100%  MeOH. 

A  second  foreign  propellant  (#2)  provided  a  greater  challenge.  After 
some  work  it  was  found  to  contain  DEGDN,  raethylcentralite*  2,6-DNT,  2,4-DNT, 
2,4,6-TNT,  nitroglycerine,  NDPA,  and  DPA.  An  internal  standard  of 
01 propylphthalate  was  selected.  A  solvent  composition  of  60%  MeOH-40%  THF  was 
chosen.  Figure  10  shows  the  resulting  chromatogram. 

A  third  foreign  propellant  (#3)  turned  out  to  be  a  single  base  nitro¬ 
cellulose  propellant  stabilized  with  diphenylaraine  (Figure  11).  The 
interesting  point  is  that  the  propellant  appears  to  be  contaminated  with  trace 
amounts  of  HMX,  RDX,  MNA,  TMETN,  4-NDPA,  and  2-NDPA.  NDPA  was  also  found  but 
It  is  probably  a  degradation  product  from  the  DPA.  Blanks  were  prepared  and 
the  contamination  was  only  found  in  the  propellant  extract.  It  was  concluded 
that  the  lannuf acturing  facility  was  also  engaged  in  production  of  other  more 
sophisticated  propellants. 

VI i.  COMPARISON  OF  TECHNIQUES 

Since  the  analysis  of  propellants  has  a  greater  confidence  level  when 
similar  values  are  obtained  by  two  different  procedures,  GC  analyses  of 
several  propellants  were  obtained  for  comparison  with  the  LC  d3ta. 
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N-5  samples  aged  In  tropic,  arctic,  or  desert  conditions  were  examined 
by  capillary  GC  for  2-NDPA  to  determine  their  remaining  shelf  life.  LC  analy¬ 
sis  of  the  samples  included  nitroglycerine,  diethylphthalate,  2,2-Dinitro  DPA, 
and  2,4-Dinitro  DPA  as  well  as  2-NDPA.  The  results  indicate  good  agreement 
for  the  two  techniques  (Table  9).  As  seen  in  the  table  and  chromatrgrams 
(Figures  12,  and  13),  formation  of  the  2-NDPA  degradation  products  proceeded 
further  , in  the  hotter  temperature  conditions  and  agrees  with  the  decrease  of 
2-NDPA.  A  potential  method  for  an  internal  check  is  to  compare  the  total 
amount  of  2-NDPA  and  its  products  with  the  initial  amount  of  2-NDPA  in  the 
propellant.  The  mlllimolar  ratios  of  2,4-Dinitro  DPA  and  2,2-Dinitro  DPA 
are  Identical  at  2.62  for  the  desert  and  tropic  samples.  The  millimolar 
amounts  of  2-NDPA  and  its  products  are  the  same  for  the  arctic  and  tropic  at 
8.37  and  8.43,  respectively.  The  desert  is  out  of  line  with  9.29  millimoles 
per  100  pg  of  propellant.  A  time  zero  analysis  for  this  lot  is  not  available; 
therefore,  it  is  difficult  to  determine  whether  the  discrepancy  is  due  to  a 
greater  amount  in  the  initial  propellant  or  an  error  in  the  analysis.  As 
pointed  out  above,  the  molar  ratios  of  dinitro  products  are  the  same  for 
desert  or  tropic  and  the  total  amount  of  2-NDPA  and  dinitro  products  is  the 
same  for  arctic  or  tropic.  Therefore,  the  most  likely  explanation  is  a  dif¬ 
ferent  amount  of  initial  2-NDPA. 

CHAPARRAL  propellant  has  been  examined  by  both  GC  and  LC  techniques  in 
RD&E  Center  laboratory.  Several  aged  motors  of  different  production  lots  were 
anatysed  to  determine  residual  levels  of  MNA,  2-NDPA,  and  4-NDPA. 

Additionally,  the  levels  of  nitroglycerine,  BTTN  and  nitroso  MNA  were  deter¬ 
mined.  Table  10  presents  a  summary  of  the  results.  Generally,  the  two  proce¬ 
dures  agree  within  a  few  percent  of  each  other  for  MNA  and  n-MNA,  but  there 
are  a  few  cases  where  as  much  as  9%  disagreement  was  observed.  The  cause  of 
this  is  not  immediately  clear.  The  total  molar  amounts  of  MNA  and  n-MNA  as 
determined  by  the  two  techniques  do  not  show  a  pattern  other  than  that  when 
aged  hot  the  LC  sees  less  material  in  the  propellant  than  the  GC  and  when  aged 
at  ambient  the  LC  see  more  material  than  the  GC.  This  may  be  caused  by  the 
incomplete  separation  of  MNA  and  HMX.  Further  study  will  be  needed  to  answer 
this  question. 

The  SPRINT  propellant  is  a  composite  modified  double-base  propellant. 
Analysis  was  performed  for  nitroglycerine,  2-NDPA,  triacetin,  resorcinol,  and 
2-nitroresorcinol .  The  latter  was  not  separated  by  GC  and  the  different 
response  factors  may  be  responsible  for  the  different  values  in  LC  and  GC. 

The  triacetin  does  not  absorb  strongly  in  the  UV  and  while  resolved  fairly 
well  small  fluctuations  in  the  base  line  or  near  by  peaks  could  significantly 
alter  the  calculated  value  (Figure  14).  Generally,  however,  there  is  good 
agreement  between  the  two  techniques  (Table  11). 

Two  foreign  propellants  were  analysed  by  HPLC  and  capillary  GC.  The 
analyses  are  given  in  Tables  12  and  13.  As  mentioned  earlier,  GC  is  not  able 
to  resolve  some  nit roderivatives  from  their  parent  compounds  where  LC  can. 

Here  NDPA  and  DPA  present  such  a  situation.  In  Foreign  Propellant  #1,  GC 
found  0.15%' DPA  while  LC  found  no  DPA  and  0.40%  nitroso  DPA.  In  Foreign 
Propellant  #2,  GC  found  0.47%  DPA  while  LC  found  0.36%  DPA  and  0.25%  NDPA 
which  corresponds,  after  correcting  for  molecular  weight  differences,  to  0.57% 
DPA  in  the  initial  propellant.  Again,  there  is  otherwise  generally  good 
agreeme  ,t  between  the  LC  and  GC  techniques. 


M-36  propellant  points  out  one  of  the  important  limitations  of  LC.  Even 
if  it  nan  be  separated,  measurement  is  impossible  if  it  is  not  visible. 
Currently,  liquid  chromatography  has  only  a  few  detection  methods  available. 

UV  is  the  most  popular,  but  many  compounds  of  interest  do  not  absorb  well  so 
they  are  detected  poorly.  Dipropyladipate  is  one  such  compound  and  as 
discussed  above  triacetin  is  another.  The  chromatogram  in  Figure  15  indicates 
that  dipropyladipate  can  indeed  be  well  resolved  from  the  other  ingredients  in 
M-36  propellant.  Table  14  indicates  that  while  several  GC  techniques  using 
flame  ionization  detectors  obtained  an  average  value  of  3.02%  di-n- 
propyladipate  with  0.5%  relative  standard  deviation,  LC  came  up  with  4.04%. 
This  is  more  than  a  33%  difference  between  techniques.  It  is  hoped  that  newer 
LC  detectors  such  as  electrochemical  or  micro  LC/Flame  detectors  can  solve 
this  problem. 

VIII.  UNIVERSAL  PROCEDURE 

No  truly  universal  procedure  exists.  In  a  production  facility  such  as 
Radford  AAP  several  different  formulations  are  being  manufactured 
simultaneously.  These  formulations  are  chemically  dif^rent  only  in  the 
amount  of  a  limited  number  of  ingredients.  Typically  s<  eral  solvent  com¬ 
positions  would  be  required  for  HPLC  analysis  of  these  fo  m  lations.  It  is 
very  probable  that  with  the  above  described  procedure  one  solvent  composition 
could  be  quickly  selected  that  would  be  able  to  separate  all, the  materials  in 
all  the  forrtulations.  Thus,  samples  could  be  analysed  as  received  or  in  a 
random  order  rather  than  holding  a  sample  until  a  group  using  the  same  solvent 
is  assembled.  As  an  example  of  this,  a  solvent  was  selected  that  could 
separate  20  of  the  most  commonly  used  propellant  ingredients  (see  Figure  7). 

It  is  clear  that  only  one  set  of  conditions  is  needed  to  quantify  any  com¬ 
bination  of  these  materials. 

IX.  CONCLUSIONS 

The  procedure  described  here  seems  to  have  almost  universal  application 
in  liquid  chromatography.  Once  the  retention  times  for  a  list  of  ingredients 
is  developed  no  further  experimentation  is  required.  Whan  a  new  propellant 
composition  is  encountered,  the  solvent  composition  necessary  for  optimum 
separation  of  the  ingredients  can  be  quickly  predicted  with  the  aid  of  a 
computer. 

The  best  separation  achieved  to  date  was  with  a  5  micrometer  spherical 
particle  coated  with  a  "monolayer"  of  C-18.  Other  columns  of  this  type  need 
to be  evaluated  to  determine  column  to  column  reproducibility.  The  computer 
program  can  evaluate  up  to  35  double  base  ingredients  but  with  the  columns 
evaluated  to  date  a  practical  limit  is  20  ingredients  per  sample.  Th*s  may  be 
improved  by  using  3  micrometer  columns.  Employing  a  second  selective  detector 
such  as  a  fluorescence  or  an  electrochemical  detector  should  offer  the  better 
solution  however. 

In  a  production  situation  this  procedure  can  quickly  predict  a  single 
solvent  composition  that  can  separate  up  to  20  ingredients.  This  should  per¬ 
mit  the  analysis  of  all  the  different  formulations  manufactured  at  a  facility 
with  one  set  of  conditions  resulting  in  a  more  complete  analysis  and  more 
rapid  turn  around  time  in  situations  presently  requiring  several  sets  of  con¬ 
ditions  or  gradient  elutions. 


It  has  been  demonstrated  that  when  a  new  formulation  of  known  composition 
is  encountered,  an  optimum  solvent  composition  can  be  predicted  without 
experimentation  once  a  column  has  beer,  characterized.  It  is  recommended  that 
when  a  column  is  replaced  with  a  similar  column  fiom  the  same  manufacturer, 
the  column  should  be  tested  at  the  optimum  solvent  composition  for  the  first 
column,  tf  no  significant  differences  are  observed  then  the  same  data  base 
can  be  used  as  before.  Columns  from  different  manufacturers  will  probably 
have  different  characteristics  and  should  be  completely  characterized 
routinely. 

Clearly,  a  procedure  exists  that  can  quickly  obtain  good  separations  of 
the  extractable  Ingredients  from  a  double-base  or  composite  modified  double¬ 
base  propellant.  LC  is  capable  of  resolving  some  materials  and  their  degrada¬ 
tion  products  that  GC  can  not.  Comparison  of  GC  and  LC  analyses  of  several 
propellants  Indicates  that  the  two  techniques  give  approximately  the  same 
results  for  those  ingredients  that  they  can  both  detect.  LC  is  superior  to  GC 
in  its  ability  to  separate  certain  materials  and  to  detect  many  nonvolatile 
materials,  such  as  HMX.  GC  is  superior  to  LC  for  detection  of  some  materials. 
LC  and  GC  techniques  thus  appear  to  complement  each  others  weaknesses  and 
strengthen  the  validity  of  the  propellant  analysis.  Therefore,  it  is  recom¬ 
mended  both  techniques  be  used  for  critical  determinations.  Otherwise,  the 
analyst  should  determine  what  he  wants  from  the  analysis  and  select  the 
appropriate  technique. 
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Table 

1.  HPLC  Optimization. 

:• 

r 

INPUT  DATA 

i 

COLUMN 

:  SUPLECO  5u  C-18  SAMPLE 

GENERAL 

>: 

DATA  IS  CORRECTED  TO 

1.00  ml/min 

PEAK  WIDTHS 

IN  SECONDS 

IN  PARENTHESES 

ACNJ 

THF& 

ACN4 

ME0H4 

*■ 

MEOH 

ACN 

THF 

MEOH 

MEOH 

THF  ' 

ACNSTHF 

1 

2 

3 

4 

5 

6 

7 

i 

Solvent  Front 

1.78 

1.79 

1.79 

1.73 

1.79 

1.79 

1.79 

■ 

i 

Resorcins! 

3.17 

3.08 

3.88 

3.24 

3.51 

3.45 

3.41  ' 

(14) 

(12) 

(12) 

(12) 

(12) 

(12) 

(12) 

2 

HMX 

3.17 

5.14 

5.61 

4.25 

5.76 

4.68 

4.72 

(14) 

(17) 

(14) 

(13) 

(14) 

(13) 

(13) 

.* 

3 

RDX 

4.07 

5.14 

6.19 

4.95 

7.28 

5.19 

5.63 

(25) 

(23) 

(17) 

(13) 

(17) 

(14) 

(14) 

-■ 

4 

DEGDN 

4.64 

6.52 

7.81 

5.93 

7.63 

6.16 

6.26 

li 

(18) 

(33) 

(18) 

(14) 

(18) 

(15) 

(15) 

| 

5 

MNA 

5.14 

5.82 

5.83 

5.95 

6.26 

5.46 

5.75 

(19) 

(18) 

(15) 

(15) 

(15) 

(14) 

(14) 

6 

Nitroglycerine 

6.49 

9.08 

16.03 

8.50 

16.07 

10.15 

10.73 

„* 

(12) 

(26) 

(44) 

(20) 

(44) 

(22) 

(22) 

7 

Acaradite  I 

5.96 

4.49 

3.65 

■6.02 

4.25 

3.90 

14.53 

(24) 

(21) 

(12) 

(15) 

(14) 

(12) 

(44) 

8 

2,4,6  TNT 

6.91 

9.25 

24.06 

9.18 

15.59 

10.80 

11.24 

(24) 

(26) 

(64) 

(20) 

(44) 

(28) 

(39) 

r 

9 

BTTN 

6.95 

10.21 

16.17 

9.66 

17.23 

10.69 

11.97 

(25) 

(28) 

(44) 

(22) 

(50) 

(28) 

(39) 

10 

2.6  ONT 

7.65 

8.70 

11.63 

9.51 

11.83 

8.46 

9.47 

(35) 

(23) 

(39) 

(21) 

(39) 

(20) 

(21) 

- 

11 

Acaradite  II 

6.86 

5.21 

4.00 

7.21 

5.65 

4.43 

5.07 

(27) 

(20) 

(13) 

(18) 

(14) 

(14) 

(14) 

12 

TMETN 

9.28 

12.81 

25.44 

12.82 

24.49 

14.05 

16.01 

*  “ 

(46) 

(34) 

(65) 

(38) 

(64) 

(39) 

(40) 

13 

NOPA 

19.96 

18.61 

17.42 

24.62 

26.50 

16.95 

23.04 

(49) 

(40) 

(44) 

(64) 

.  (64) 

(40) 

(58) 

14 

4  NOPA 

17.86 

15.07 

14.84 

21.09 

22.06 

13.72 

'  18.81 

V 

(44) 

(40) 

(40) 

(53) 

(55) 

(39) 

(50) 

5 

15 

Meth  Central ite 

14.64 

9.39 

6.51 

15.32 

8.04 

7.01 

8.69 

9 

(52) 

(36) 

(17) 

(40) 

(20) 

(17)  ■ 

(22) 

9 

16 

DPA 

36.63 

19.36 

11.13 

38.97 

16.29 

12.35 

18.17 

fi 

(55) 

(44) 

(39) 

(60) 

(38) 

(36) 

(40) 

i 

17 

2  NOPA 

35.06 

24.40 

18.75 

37.27 

30.74 

19)32 

26.91 

L 

(104) 

(64) 

(40) 

(60) 

(55) 

(52) 

(68) 

L 

18 

Ethyl  Centralite 

14.89 

15.54 

6.48 

19.39 

14.59 

11.55 

14.20 

(117) 

(79) 

(15) 

(55) 

(40) 

(39) 

(40) 

19 

Triacetin 

3.88 

4.17 

3.61 

4.54 

4.05 

3.73 

3.90 

1- 

(19) 

(27) 

(49) 

(15) 

(40) 

(29) 

(33) 

9 

20 

Tetryl 

4.28 

5.24 

6,13 

5.29 

7.99 

5.24 

5.95 

S 

(21) 

(22) 

(26) 

(83) 

(84.) 

(23) 

(26) 

2 

21 

TEGON 

5.02 

7.14 

7.44 

6.62 

7.64 

5.99 

6.40 

(31) 

(21) 

(23) 

(29) 

(25)  . 

(19) 

(21) 

22 

Ni  troso-MtlA 

6.11 

7.23 

8.31 

7.62 

8.64 

6.59 

7.26 

• 

(27) 

(20) 

(27) 

(21) 

(26) 

(20) 

(23) 

23 

2  N-Resorcinal 

3.83 

4.20 

4.89 

4.35 

5.17 

4.45 

4.63 

(2D 

(14) 

(16) 

(15) 

(19) 

(15) 

(16) 

24 

2,4  ON  Resorcin 

1.89 

1 .84 

1.84 

2.00 

2.19 

1.89 

2.03 

(16) 

(16) 

(20) 

05) 

(24) 

(16) 

(20) 

25 

2,4  ONT 

8.07 

8.68 

11.61 

9.68 

13.20 

8.57 

10.04 

(33) 

(21) 

(32) 

(26) 

(36) 

(24) 

(28) 

26 

2,4,4  TN  DPA 

13.78 

15.39 

10.02 

20.82 

29.28 

16.43 

20.13 

(56) 

(42) 

(80) 

(60) 

(75) 

(54) 

(62) 

# 

27 

2,4  ON  DPA 

17.21 

16.53 

19.20 

22.75 

23.58 

14.92 

19.14 

‘ 

(53) 

(44) 

(51) 

(65) 

(70) 

(68) 

(92) 

28 

2,2  ON  0PA 

21.60 

18.03 

12.66 

26.46 

24.19 

14.08 

19.82 

(83) 

(27) 

(37) 

(73) 

(71) 

(41) 

(57) 

29 

Methylphthalate 

5.44 

5.85 

4.62 

6.64 

5.33 

4.84 

6.22 

,  - 

. 

(25) 

(17) 

(23) 

(21) 

(19) 

(24) 

(18) 

30 

Ethyl  Phthalate 

9.72 

9.82 

6.85 

12.41 

8.21 

1  7.29 

8.35 

(42) 

(27) 

(22) 

(37) 

(94) 

(33) 

(115) 

L 

31 

Propylphthalate 

24.38 

21.52 

12.36 

32.20 

17.49 

13.82 

18  07 

(93) 

(58) 

(36) 

(91) 

(56) 

(41) 

(53) 

■ 
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Table  2.  HPLC  Optimization. 

9:00  AM  TUE.,  14  AUGUST  1984 
DATA  ANALYSIS 

COLUMN:  SUPELCO  5u  C-18  SAMPLE:  N-5  PROPELLANT 

THE  BEST  COMPOSITION  IS: 

MEOH  ACN  THF 

82%  18%  0% 

THE  WORST  RESOLUTION  AT  THIS  POINT  IS  2.855 

PREDICTED  ORDER  OF  ELUTION 

RELATIVE  RETENTION  TIME 

INGREDIENT  1.00  ml/mln  1.50  ml/min 


Nitrogylcerlne 

7.38 

4.92 

Ethyl  Phthalate 

11.30 

7.54 

2,4,4  TN  DPA 

17.76 

11.84 

2,4  DN  DPA 

20.57 

13.71 

2,2  DN  DPA 

24.89 

16.59 

Propylphthalate 

29.34 

19.56 

2  NDPA 

37.61 

25.07 

■ci  «r* 


Table  3.  Optimization  of  5  Micron  Waters 
C-18  Radial  PAC. 
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MEOH 

ACM 

THF 

MEOH/'ACM 

MEOH/THF 

ACM.' THF 

IIEOH.-ACH/'THF 

TRIflCET  III 

1.60 

3.21 

8.60 

1 . 62 

7.  03 

1 . 99 

4.69 

RESORCINOL 

1.81 

3.66 

6.90 

1.91 

7.34 

2.10 

5.05 

MMX 

1 . 95 

3.32 

3.77 

2.82 

4.34 

3.  17 

3.  34 

RDX 

2. 86 

3.41 

4.33 

3.31 

6.13 

3.  70 

4 . 25 

DEG  DM 

3.  40 

4.51 

3.96 

,  4.13 

6.  19 

4.36 

4. 77 

HNR 

4.  12 

4.23 

4.30 

4.31 

5.39 

4.15 

4.66 

MC 

4. 33 

6.36 

14.45 

5.69 

12.66 

8.  12 

9.05 

ACAPADI TE 

I  4.74 

3.20 

2.33 

4.40 

3.15 

2.70 

3.  2 7* 

TNT 

4.63 

6.42 

21.72 

6.16 

13.30 

6.  78 

9.00 

DTTM 

4.88 

7.  12 

13.40 

6.65 

13.66 

8.33 

9.  77 

2,6  put 

3.93 

6.26 

9.70 

6.  58 

?.  96 

7.  20 

8.26 

ACAF:AD  1  TE 

I  I  6.19 

4.11 

2.66 

3.35 

3.62 

3.21 

3.65 

TMETH 

6.43 

9.00 

22.05 

8.68 

22.09 

11.41' 

13.28 

2, 4  DMT 

6.  93 

6.36 

9.63 

7.03 

10.50 

7.  12 

6.41 

EC 

12.13 

11.21 

10.10 

13.  30 

12.98 

10.51 

13.45 

MC 

13.14 

9.  C9 

? .  93 

12.02 

6.60 

5.63 

7.67 

A  NDFA 

13.70 

11.22 

12.43 

15.81 

22.00 

11.63 

17.28 

IIPFA 

13.  96 

12.96 

15.  95 

13.78 

23.30 

15.50 

1  21.60 

£  MOPA 

29.  60 

17.81 

17.20 

26.70. 

26.45 

17.52 

27.  15 

DPR 

31.30 

16.40 

10.13 

27.  88 

14.85 

11.40 

17.98 

Table  4.  Optimization  of  10  Micron  Waters 
C-18  Radial  PAC. 


MEOH 

AC  li 

THF 

MEOH. ACN 

MEOH.'  THF 

AC  II.' THF 

ME  OH 'ACM.' 

TF 1ACETIM 

.93 

1.16 

1 . 59 

1 . 38 

1.00 

1 . 38 

1 . 50 

PE  sore  I  MOL 

.66 

.  30 

1 . 05 

.  90 

1.03 

.93 

.  93 

HM}{ 

.  93 

1  .  40 

1 . 57 

1  .  40 

1.87 

1 . 37 

1.51 

FOX 

1.40 

1 . 33 

1.82 

1 . 69 

2.49 

1 . 57 

1  .  66 

DECDM 

1 . 75 

2.01 

2 . 

2.13 

2.71 

1 . 99 

2. 30 

Mil  A 

2.  10 

1 . 63 

1 . 82 

2.  34 

2.30 

1 . 80 

2.  10 

MC 

1 

2 .  £ 

1.17 

3.06 

5.  92 

3*38 

5  f  Ci* 

acafamte  J 

1.45 

1.01 

2.30 

1 . 37 

1.16 

1.45 

TNT 

2.54 

2. 94 

9.  D? 

3.  38 

►  .  0* 

3 ,  r 

4 . 06 

6T  Til 

2. 57 

3.  10 

3.1  1 

3.71 

6.61 

2.  f  3 

4.43 

2,6  MIT 

3.17 

2.80 

4. 09 

3.53 

4.46 

£.91 

3.4? 

ACAFAD1TE  11 

7. 

1 . 89 

1.15 

3.  02 

1.57 

1 :37 

1.71 

TME  tm 

3.54 

4  .  0* 

1  7. 58 

5.02 

?  .  t  6 

5.00 

6.16 

2,4  MIT 

6.93 

€  .  :*«? 

9 .  ► 

r .  0  ;• 

1  0.  so 

7.  12 

6.4  1 

EC 

19.36 

r .  £?• 

4 . 08 

16.51 

6.51 

4 . 30 

7.15 

MC 

7.65 

3.68 

2.  1  1 

e.  f  t 

2.96 

»■  «T 

a.  •  »■  w> 

3. 24 

4  NO  FA 

7.59 

5.  12 

5.  74 

8.97, 

9.17 

5.  1  7 

7  .  4  6 

lipFA 

7.09 

3.13 

6.4  1 

7. 68 

5.78 

4.17 

3. 34 

2  NOP  A 

18.17 

8.  24 

7.11 

15.37 

12.  91 

7.18 

11.22 

IiF  A 

7.  66 

6.4  1 

9.  00 

10.34 

6.14 

6. 75 

11 
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Table  5.  Optimization  of  5  Micron  Perkin-Elmer 
HS-5  C-18  Column. 


MEOH 

RCM 

THE 

MEOH'-RCM 

MEOH/IHF 

RCM'THF 

MEDH^RCN^THF 

TF:1RCET1M 

2.27 

2.31 

2.96 

2.43 

2.99 

2.75 

2.68 

RESORC1M0L 

2.27 

2.31 

2.96 

2.43 

2.99 

*>.  75 

2.66 

MMX 

2.54 

4.21 

5.63 

3.58 

5.22 

4.28 

4.09 

RPX 

3.37 

4.51 

6.  34 

4.15  . 

6.  75 

4.78 

4.84 

PEGDH 

3.85 

5.47 

8.84 

4.89 

6.84 

5.63 

5.42 

mmr 

4.53 

5.07 

5.  84 

5.  14 

5.91 

4.96 

5.09 

HC 

5.  14 

8.  18 

17. 00 

6.91 

14.51 

9.92 

9.78 

RCRRRD  1  TE  1 

4.54 

3.48 

3.  13 

4.69 

3.  72 

3.34 

3.  69 

TUT 

5.42 

8.27 

29.00 

7.24 

14.44  . 

10.42 

9.43 

PTTM 

5.56 

0.22 

18.75 

8.33 

ie.es 

10.43 

10.80 

2,6  DMT 

6.27 

7.69 

12.26 

7.65 

10.69 

6.03 

8.36 

RCRRflPlTE  11 

6.12 

4.38 

3.  57 

5.  79 

4:i9 

3.  92 

4.34 

TMETM 

7.46 

11.02 

30.00 

10.77 

23.21 

13.84 

14.40 

2,4  DMT 

6.58 

7.63 

11.88 

7.51 

10.  79 

8.  1  1 

e.  29 

EC 

12.  70 

12.94 

11.86 

14.83 

13.59 

11.14 

12.56 

ric 

13.31 

8.  18 

6.40 

11.93 

7.69 

6  3 

7.89 

4  MPRR 

14.55 

12.75 

17.  34 

17.C6 

22.  20 

13.  33 

16. 65 

MDFR 

14.94 

14.93 

19.26 

17.56 

38.23 

15.91 

19.35 

2  MDFR 

27.86 

18.96 

20.99 

26.54 

28.  22 

16.05 

23.17 

DPR 

34.. 73 

15.  75 

11.86 

29.06 

16.02 

11.96 

16.30 

Table  6.  Optimization  of  5  Micron  Supelcosil 
LC-18  Column. 


MEOH 

RCM 

THF 

MEC'K-'RCM 

IIEOH^THF 

RCM.'  T  HF 

ME  OH/'RCM''  THF 

TRIRCETIM 

3.  1  7 

3.  OS 

3.88 

3. 24 

3.51 

3.45 

3.41 

RESORCINOL 

3.  17 

3.08 

3.68 

3.24 

3.51 

3.45 

3.41 

HMX 

3.  1  7 

5.  14 

5.61 

4.25 

5.  76 

4.68 

4. 72 

RDX 

4.07 

5.14 

6.  19 

4.95 

7.28 

5.  19 

5.63 

PE  COM 

4.64 

6.52 

7.61 

5.93 

7.63 

6.  16 

6.26 

MMR 

5.  14 

5.  62 

5.  e? 

5.95 

6.26 

5.46 

5.  75 

MG 

6.49 

6.08 

16.03 

6.50 

16.07 

10.  15 

10.  73 

rcrrrdite  1 

5.  96 

4. 49 

3.65 

6.02 

4.25 

3.90 

14. 53 

TMT 

6.91 

9.  25 

24.06 

9.  16 

15.59 

10.60 

11.24 

PTTM 

6.95 

10.21 

16.17 

9.66 

17.22 

10.69 

11.97 

2,6  PMT 

7.65 

6.  70 

11.63 

9.51 

11.63 

6.46 

9.47 

RCRFRDITE  11 

6.  86 

5.21 

4.06 

7.21 

4.65 

4.43 

5. 07 

TMETM 

9. 28 

12.81 

25.4  4 

12.62 

24.49 

14.65 

16.01 

2,4  DMT 

8.03 

8.71 

11.45 

9.47 

•2.08 

6.75 

9.  38 

EC 

14.89 

15.54 

6.43 

19.39 

14.5? 

I  1 .55 

14.20 

MC 

14.64 

9.  39 

6.51 

15.  32 

8.64 

7.01 

6.  69 

4  MDFR 

17. 86 

15.67 

14.64 

21.09 

22.06 

13.72 

18.61 

MPF  R 

19.  96 

13.61 

17.42 

24.62 

26.  50 

16.95 

23.04 

2  MDFR 

35.06 

24. 40 

18.  75 

37.27 

30.  74 

19.  32 

26.91 

DPR 

36.63 

19.  36 

11.13 

36. 97 

16.29 

12.35 

18.  17 
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table  7.  Optimization  of  5  Micron  Dupont  Zorbax  ODS  Column 


HE  OH 

ACN 

THE 

NEOH/ftCN 

HEDH/THF 

ACIWHF 

tlEOH'ACN/'THF 

TR1ACET IN 

3.06 

2.08 

3.53 

3.  30 

3.50 

3.42 

3. 49 

RESORCINOL 

3.06 

2.08 

3. 53 

3.  30 

3.50 

3.  42 

3.49 

HMX 

3.23 

5.2  2' 

5.48 

4 . 55 

5. 88 

4.88 

5.  1  7 

RDX 

s.  re 

5.56 

6.23 

5.  64 

8.69 

5.53 

6.39 

ItEGDN 

6.19 

7.11 

8.01 

7.33 

8.74 

7.  07 

7.49 

TINA 

7.33 

6.56 

6.09 

7.74 

7.47 

6.31 

6.89 

NG 

F.  49 

10.47 

20.13 

1  0.  38 

22.24 

12.43 

13.82 

AC  ARAM  TE 

I  9.07 

5.56 

3.  52 

7.07 

4.60 

4.06 

4  .  ee 

TNT 

9.67 

10.95 

29.00 

11.64 

21.79 

1  3 .  F  3 

14.75 

t'TTN 

9.  74 

11.57 

19.77 

12.04 

23. 28 

13.45 

15.43 

2,6  It  NT 

11.03 

10.53 

13.21 

12.14 

16.39 

10.  05 

12.35 

ACARAI'l  TE 

1111.53 

7.  1  1 

3.96  8 

10.03 

5.19 

4 . 72 

5.65 

THE  Til 

13.31 

1  4  . 

30 . 50 

16.90 

32.26 

IF.  30 

21.60 

2,4  I' NT 

14.15 

1  0.  SI 

1  3.  E  2 

13.12 

17.47 

10.84 

12.62 

EC 

27.  12 

20.4  6 

12.32 

26. 95 

20.  71 

14.31- 

19.17 

nc 

31.35 

16.11 

7. 03 

25:6  2 

10.  28 

8.25 

10.95 

4  NI'FR 

If.  48 

18.76 

17.76 

30. 00 

31 . 00 

16.27 

25. 32 

NO  FA 

31.92 

24  .ft? 

re*.  f-c' 

33.  35 

39.  1  7 

21 . 9F 

31.97 

2  NltFft 

7  ? .  S3 

34040 

24.  ^E 

59.  77 

4F.  35 

2  7.47 

40.31 

It  F  A 

70.  04 

34. 61 

m 

f ...» 

62. 94 

22.35  . 

15.31 

24.15 

Table  8,  Column  Specifications. 


Source 

nimens ions 

Packing 

Partical 

(MM) 

Plates 


Waters 

100 

X 

5 

u-Bondapac  C-18 

10 

Micron 

Spherical 

>5,000 

Waters 

100 

X 

5 

u-Bondapac  C-18 

5 

Micron 

Spherical 

>9,000 

Supelco 

250 

X 

4.6 

Supelcosil  EC-18 

5 

Micron 

Spherical 

20,614 

E  Merck 

2  50 

X 

4.6 

LiChrosorb  RP-18 

5 

Micron 

Irregular 

15,516 

Dupont 

250 

X 

4.6 

Zorbax  ODS 

5 

Micron 

Spherical 

27,135 

PERKIN- 

E1..MEH 

125 

X 

4.6 

P-E  US—  5  C-18 

5 

Micron 

Spherical 

>8,000 
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Table  9 

.  N-5  Propellant  Analysis. 

NG 

DEP 

2,2-DN  DPA 

2,4-DN  DPA 

LC 

2-NDPA 

GC 

2-NDPA 

Arctic 

35.74 

10.72 

— 

— 

1.792 

1.811 

Tropic 

35.88 

10.88 

0.090 

0.236 

1.534 

1.716 

Desert 

35.91 

11.13 

0.201 

0.526 

1.388 

1.379 

Table  10. 

CHAPARRAL  Propellant  Analysis. 

Comparison  of  Capillary  GC  and  HPLC 

Lot 

Condi tion 
Aged 

Analytical 

Method 

MNA 

N-MNA 

2-NDPA 

•  4-NDPA 

NG 

BTTN 

E0001 

Hot 

GC 

0.259 

0.490 

0,244  ' 

0.197 

12.56 

12.86 

Hot 

LC 

0.263 

0.334 

0.245 

0.209 

13.43 

12.23 

E0003 

Ambient 

GC 

0.564 

0.043 

.0.260 

0.229 

13.18 

13.55 

Ambient 

LC 

0.609 

0.033 

0.257 

0.230 

■  13.14 

13.45 

E0003 

Hot 

GC 

0.346 

0.380 

0.255 

.  0.223 

13.04  ' 

13.26 

Hot 

LC 

0.366 

0.281 

0.244 

0.216 

13.51 

13.62 

E0005 

Ambient 

GC 

0.588 

0.048 

0.260 

0.248 

12.76 

12.79 

Ambient 

LC 

0.622 

0.029 

0.256 

0.256 

12.84 

12.74 

E0005 

Hot 

GC  . 

0.112 

0.727 

0.241 

0.211 

12.38 

14.80 

Hot 

LC 

0.107 

0.624 

0.235 

0.204 

13.00 

14.99 

E0006 

Ambient 

GC 

0.569 

0.042 

0.251 

0.222 

12.81 

14.57 

Ambient 

LC 

0.623 

0.022 

0.239 

0.220 

12.75 

14.45 

E0007 

Ambi ent 

GC  • 

0.597 

0.C39 

0.260 

0.250 

12.95 

12.97 

Ambi  ent 

LC 

0.637 

0.024 

0.247 

0.249 

12.80 

12.78 

E0008 

Ambient 

GC 

0.626 

0.023 

0.261 

0.254 

13.26 

12.64 

Ambient 

LC 

0.682 

0.022 

0.254 

0.250 

12.94 

12.46 

E0008 

Hot 

GC 

0.205 

0.644 

0.340 

0.218 

13.03 

13.65 

Hot 

LC 

0.218 

0.530 

0.323 

0.210 

13.4V 

13.73 

GCU 

Ambient 

GC 

0.615 

0.115 

0.256 

0.261 

13.28 

14.13 

Ambient 

LC 

0.659 

0.061 

0.248 

0.239 

13.66 

14.34 
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Table  11.  Sprint  Propellant  Analysis. 


Triacetln 

Resorcinol 

2-Nit rorescinol 

NG 

2-NDPA 

LC 

6.657 

0.950 

0.019 

28.734 

0.942 

GC 

6.325 

0.901 

28.400 

0.966 

• 

Table  12. 

NG 

Foreign  Propellant 

DPA 

#1  Analysis. 

MC 

NDPA 

GC 

39.25 

0.153 

1.104 

— 

LC 

39.41 

1.345 

0.404 

Table 

13.  Foreign  Propellant  #2  Analysis. 

DEGDN 

NG 

2,6  DNT  +  2,4  DNT  2,4,6  TNT 

MC 

DPA 

NDPA 

LC 

14.22 

17.75 

(0.905)  0.375 

2.516 

0.361 

0.249 

GC 

14.64 

17.41 

(0.918)  0.356 

2.550 

0.474 

— 

Table  14.  M-36 

Propellant  Analysis. 

NG 

DNDPA 

2-NDPN 

GC 

Splitless 

39.151 

3.017 

1.950 

GC 

On-Coluran 

39.134 

3.046 

1.961 

GC 

Manual/Split 

39.117 

3.011 

1.969 

GC 

Auto/Split 

39.182 

3.002 

1.962 

GC 

Auto/Split 

39.012 

3.016 

1.954 

LC 

38.626 

4.040 

1.777 

/ 
/  ■ 
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Figure  2.  Standard  sample. 


10  MICRON  WATERS  BEST 


Figure  4.  Waters  5  micrometer  radial  FAC  column. 


Figure  5^  Waters  10  micrometer  radial  PAC  column. 
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Figure  7.  Suplecosil  LC-18  column. 


igure  9.  Waters  10  micrometer  radial  Pac  column. 
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ellant  #3 
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Figure  14.  SPRP^  propellant  standard. 


M-36  PROPELLANT 


a 


LIST  OF  ABBREVIATIONS 


HPLC 

GC 

MeOH 

ACN 

THF 

DPA 

NDPA 

2-NDPA 

4-NDPA 

2,2-DNDPA 

2.4  "-DNDPA 

2.4.4  ''-TNDPA 
2,4-DNT 

2.6- DNT 

2 .4 .6- TNT 
BTTN 
DEGDN 
DNDPA 
MNA 
NMNA 

NG 

TIMETN 


High  Pressure  Liquid  Chromatograph}' 

Gas  Chromatography 

60%  Methanol  WiCer  Mixture 

50%  Acetonitrile  Water  Mixture 

42%  Tetrahydrofuran  Water  Mixture 

Diphenylamine 

Nitrosodiphenylamine 

2-Nitrodiphenylamine 

4-Nitrodiphenylamlne 

2,2-Dinitrodiphenylamine 

2 .4  "-Dinitrodipher.y!  amine 

2.4.4  "-Trinitrodiphenylamine 

2 . 4- Dini t rotoluene 

2.6- Dinitrotoluene 

2.4.6- Trinitrotoluene 

1 .2.4- Butanetrioltrinitrate 
Diethyleneglycoldinitrate 
Di-n-propyladipate 

N-M ^thylanaline 

N-Nl  troso-n-methylanaline 

Nitroglycerine 

2-hydroxymethyl-2-methyl-l , 3-propanediol trinit rate 
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The  analog  out  put  from  the  detector  was  digitized  and  sent  to  a 
Hewlett-Packard  Laboratory  Automation  System.  The  data  was  then  integrated  in 
the  normal  manner.  A  post  run  basic  language  program  "PEAKWT"  listed  below 
was  then  automatically  run.  This  program  examines  the  data  file  for  the  last 
chromatograph.  For  each  peak  integrated,  the  p-ogram  finds  the  maximum  ampli¬ 
tude.  It  then  samples  the  data  points  before  and  after  the  peak  until  it  has 
a  corrected  value  of  10%  of  the  maximum.  The  retention  time  of  these  points 
are  then  subtracted  and  reported  as  the  width  of  that  peak. 

The  10%  was  arbitarlly  selected  because  it  is  more  reproducible  and  does 
not  confuse  a  tailing  peak  with  a  very  broad  peak.  This  system  gives  more 
meaningful  results  when  all  ingredients  are  well  separated  and  of  approxima¬ 
tely  the  same  amount. 

In  Line  190  A$  will  be  returned  string  with  the  name  of  the  raw  data  file 
and  B$  the  name  of  the  processed  data  file.  In  line  240  and  250  W  is  the 
width  of  a  data  slice  and  N  is  the  number  of  peaks  integrated.  In  line  290 
T(I)  and  P(I)  are  the  retention  time  and  amount  of  pc.ak  I.  Line  320  sets  the 
pointer  to  the  first  record  containing  extended  data.  Line  350  reads  the 
height,  start,  baseline  at  start,  end,  and  baseline  at  end  of  each  peak  I. 

The  Loop  in  480  to  540  finds  the  10%  levels  at  the  front  and  end  of  each  peak 

I. 


100  REM  *»••*•»*•«*•*»»»*»**•***»****»•*******•»*****•*«*»»*»•******* 

110  REM  *  *  *  BASIC  PROGRAM  "PEAKWT"  **» 

120  REM  **»  CALCULATES  PEAK  WIDTH  AT  10*  OF  THE  PEAK  HEIGHT.  » * » 

130  REM  ***  TO  CALCULATE  AT  A  DIFFERENT  FRACTION  OF  THE  PEAK1  **» 

140  REM  «•*  HEIGHT  CHANGE  LINES  360  AND  370.  *** 

150  REM  ***  CREATED  BY  J.  CARVER  7  JUNE  1984  *** 

160  REM  ***•**»»•»••*•«**•****•*»•»**•*****«**»»♦*♦*»*«**•**»•*•*•»»« 
16  5  COM  35(101 ,C$(1Q1  ,MS[10]  ,RS[10)  ,P$(101 ,N|101  ,0(2,41 ,Z(4| 


170  DIM  NS(261  ,A(35]  ,3  (351  ,C'  35]  ,E'35]  ,H  (351  ,P(35!  ,S(35'i  ,T(  351  ,U(351 
180  FILES  *,* 

190  TELL  A S,0 $  ,C S  ,S 
200  ASSIGN  8S  ,  1 , E , F 
210  ASSIGN  CS, 2 , E  ‘ 

220  ON  END  II  THEN  1000 
2  3 J  ON  END  12  THEN  1000 

24  0  READ  12,  Is  NS ,X ,X ,N S ,X ,NS ,X ,X  ,X  ,X  ,N S ,N S ,N S , NS ,NS ,X ,X ,W 
2  50  READ  *1,1: N$ ,X,X .NS.X ,NS ,X ,X ,X,X ,N$,NS,N3,NS ,NS,X ,X ,X ,X ,X  ,M 
260  LET  N=ABS  (N) 

270  READ  41,2 
280  FOR  1*1  TO  N 

290  READ  # 1; T  1 1 )  .X,X,X,X,P[  I)  ,N$ 

300  NEXT  I 

310  LET  R4*l 
320  READ  II, ( F/2) +2 
330  READ  12,2 

FOP  1*1  TO  N 

READ  I  1:  H  ( 1 1  ,S  ( I  )  ,  X  ,B  (  I  ]  ,  E  [  I  1  ,X,U(I)  ,NS 
LET  A  (I  1= (H ( I  1* .1) +B(I 1 
LET  C(I1*(II[I]*.1)+U(I1 
LET  S[I1*S(I!»60 
LET  E(I|*E[I]*60 
LET  Tfl] «T (11*60 
LET  P.  1  *1  NT  ( S  ( 1 1  /W) 

LET  R2*INT (E (I  1 /W) 

LET  P.3  =  IMT(T(1 1 /W) 

FOR  K*R4  TO  Rl-1 
READ  42: N2 
NEXT  K 

IF  R4  >*  R 1  THEN  LET  R1-R4 
FOR  K*R1  TO  R2 
READ  42 ;N2 
IF  K>  P.3  GOTO  530 

IF  M2  >*  A ( I ]  AND  INT (S ( I  1 /W) «R1  THEN  LET  S(I]=K*W 
GOTO  540 

IF  N2  <*  CU1  AND  INT  ( E  ( 1 1  /W)  *R.  2  THEN  LET  E(I)*K*W 
NEXT  K 
LET  R4  =R2  +  1 
NEXT  I 
1000  PRI NT 
1010  PRINT 

1020  PRINT  " PPOCESSED  FILE  ";BS 
1030  PRINT 

1040  PRINT  TA  B  (  30 )  ; " 10 *  OF  HEIGHT" 

1050  PRINT  "PEAK  RT" , "AREA  %  "."PEAK  WIDTH  (SEC)" 

1060  PRINT 

1070  FOR  1*1  TO  N 

1080  LET  W=(E (I  1 -S (I  1 ) 

10  90  LET  T ( I  1 =T ( I  1 /60 

1100  PRINT  T (  1 1  ,  P (  I  1  ,W 

1110  NEXT  I 

1120  PPINT 

1130  PFINT 

1140  PRINT 

1150  END 


34  0 
350 
360 
370 
380 
390 
400 
410 
4  20 
430 
440 
450 
460 
4  70 
480 
490 
500 
510 
520 
530 
540 
550 
560 


A-3(A-4  Blank) 


APPENDIX  B 

DATA  REDUCTION  BASIC  LANGUAGE  VERSION 


* 

The  basic  language  version  of  the  data  reduction  was  designed  to  output  to 
a  Hewlett-Packard  2648  graphics  terminal  with  a  9876  A  graphics  printer 
attached  as  a  "C"  device. 

Line  145  to  155  In  program  XLDB6A  open  previously  created  data  files  to 
the  program.  These  7  files  contain  the  retention  data  from  the  7  eolvent  com¬ 
positions.  This  program  can  enter  new  data  or  change  old  data  in  these  files. 
It  can  not  create  the  files.  This  must  be  done  separately  by  the  user^  When 
the  files  are  satisfactory,  the  user  can  list  the  data  input  as  well  as  the 
resolution  for  pairs  of  peaks  in  each  solvent.  A  chromatographic  optimization 
function,  as  discussed  by  Synder  [21],  is  also  calculated.  The  program  then 
calls  XLDB6B  for  further  calculations  and  plotting  lines  of  resolution  for 
each  peak  pair  on  the  graphics  terminal. 


10  COM  M|2J,T(7,35] 


100  REM  ••***•••*••*****•***••*•****«•••»•*****»***• 
105  REM  •*»***»**MQUIO  CHROMATOGRAPHIC  OPTIMIZATION 
110  REM  *»••******•  XLDB6A  *  *  R . RADKE  /  J. CARVER  **** 
115  REM  **********«*****»»*«*****aUG/SEPT  1981****** 
120  REM  •**•**»**••****«***•«•••••*•*•*•*•*•*******• 


125  LET  Ml«7  , 

130  LET  M2-8 

135  DIM  E$|l] ,S{71 ,V(7) ,W(7,35| ,K| 7,35! ,R[7, 35) 

140  CHRS  27, ES 

145  PILES  FIRST: :19, SECOND: : 19 ,TH I RD : : 1 9 
150  TILES  POORTII:  :19, FIFTH:  :  19, SIXTH  ::  19 
155  FILES  SEVNTII :  :  19 
160  REM  GET  SOLVENTS 
165  PRINT 

170  PRINT  *  ENTER  THE  NUMBER  OF  SOLVENTS  (1-7)  <0-EXIT>" 

175  INPUT  Ml 

180  IF  Ml»0  THEN  9999 

185  IF  M1<1  THEN  170 

190  IF  Ml>7  THEN  170 

195  LET  M { 1 ) *M l 

200  REM  GET  PEAKS 

205  PRINT  *  ENTER  THE  NUMBER  OF  PEAKS  PER  SOLVENT  (2-35)  <0»EXIT>" 

210  INPUT  M2 

215  IF  M2-0  THEN  9999 

220  IF  M2< 2  THEN  205 

225  IF  M2>35  THEN  205 

230  LeT  M(2|-M2 

235  REM  LOOP  TO  GET  DATA 

240  FOR  J»1  TO  Ml 

245  PRINT  *  ENTER  THE  RETENTION  TIME  (MIN)  FOR  SOLVENT  "jJ;*  <0-EXtT>* 

250  INPUT  V(J| 

255  IF  V(Jl-0  THEN  9999 

260  LET  V (J | «V (J ) *60 

265  PRINT  "ENTER  THE  AVERAGE  BAND  WIDTH  <SEC>  FOR  ALL  PEAKS  OF  SOLVENT  *;J 

270  INPUT  A 

275  FOR  1-1  TO  M2 

280  PRINT  *  ENTER  THE  RETENTION  TIME  (MIN)  FOR  PEAK  *jlf*  <0-EXIT>" 

285  INPUT  T ( J , I ) 

290  IF  T|J,I]»0  THEN  9999 

295  LET  T(J,I )-t:J.I]*60 

300  LET  W[J,l )«A 

305  LET  K f J , I J « (T I J , I | -V | J | ) /V ( J ) 

310  NEXT , I 

315  NEXT  J 

320  GOSUB  605 

325  REM  CALCULATE  THE  RESOLUTION 

330  GOSUB  6000 

335  FOR  J*=l  TO  Ml 

340  FOR  l -1  TO  M2-1 

345  FOR  K-I+l  TO  M2 

350  LET  A“ (T | J ,K|~T( J , I ! ) / ( . 5* (W |J , K | *W (J , I j ) ) 

355  PRINT  IJjA 

360  NEXT  K 

365  NEXT  I 

370  NEXT  J 

375  REM  •*»  CALCULATE  RESPONSE  FUNCTION 

3R0  PRINT  "  ENTER  DESIRED  RESOLUTION  <0»GO  TO  ORM>" 

385  INPUT  R9 

390  IF  R9»0  THEN  CHAIN  " XLDB6B" 
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395 

IP  R9<0  TURN  380 

400 

GOSUB  6000 

405 

FOR  J-l  TO  Ml 

410 

LET  S|J)«0 

415 

FOR  l-l  TO  M2-1 

420 

FOR  K-I+l  TO  M2 

425 

REAP  #J;A 

4  30 

LET  Z-A/R9 

435 

IF  Z  <«  0  Til FN  445 

44  0 

LET  S{J)«S|J)*LN(Z) 

445 

NEXT  X 

450 

NEXT  I 

455 

NEXf  J 

4*0 

REM  ***  OUTPUT  TABLE  **• 

4*5 

GOSUB  6000 

470 

FOR  J«1  TO  Ml 

475 

OUTPVC  "Cl",E 

480 

PRINT  E$"U* 

485 

PRINT  *  RETENTION  TIME  FOR  SOLVENT  ";J;" 

IS  *;  (VUI/60)  ;*  MIN" 

490 

PRINT 

495 

PRINT 

500 

PRINT  *  PEAR*,*  RETEN  TIME*,"  BAND  WIDTH 

","  CAPACITY* 

505 

FOR  !«1  TO  M2 

510 

PRINT  TAB (3) jl jTAB ( 19) ; (T | J, I | /60) i TAB  I 34 ) ; Wl J , I | ;TAB(47) ;K|J 

515 

NEXT  I 

520 

PRINT 

525 

PRINT 

5  30 

PRINT  *  FOR  PEAKS',"  RESOLUTION" 

535 

PPINT 

540 

FOR  I»1  TO  M2-1 

545 

FOR  K» I ♦  1  TO  M2 

1 

550 

REAP  f J ; A 

555 

PRINT  TAB <  3)  ;  1 J "  “  ";K;TAB{ 17) ;A 

5*0 

NEXT  K 

56  5 

NEXT  I 

570 

PRINT 

575 

PRINT 

580 

PRINT  "  CHROMATOGRAPHIC  OPTOMIZATION  FUNCTION  (COF)  FOR  SOLVENT 

58  5 

PRINT  *  IS  ";S|J1 

590 

NEXT  J 

595 

OUTPVC  "  T  *  "  ,  E 

6  00 

CHAIN  "XLPB6B" 

*05 

REM  PPINT  INPUT 

610 

PPINT 

615 

PRINT 

620 

PRINT  TA B ( 8 ) “ SOLVENT- 1  SOLVENT- 2  SOLVENT-3 

SOLVENT-4  SOLVENT-5"; 

6  25  , 

PRINT  TAB  (58)  "SOLVENT-*  SOLVF.NT-7* 

O  Ju 

635 

PRINT  TAIM58)  * - ---" 

640 

FRINT  "TO"  ; 

64  5 

FOP  X  *1  TO  Ml 

650 

PRINT  TAB(X»10-2) ;V(X|/60; 

655 

NF  XT  X 

660 

PRINT 

665 

PRINT  "WIDTH"; 

670 

FOR  X* l  TO  Ml 

■ 

675 

PRINT  TAB (X* 10-2)  ;W(X  ,1)  ; 

680 

NEXT  X 

*85 

PPINT 

690 

fOR  Y*  1  TO  M2 
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*•**••>  PRINT  *RT-"Y; 

/no  FOR  x-1  TO  Ml 

70 S  PRINT  TAIMXMO-2)  ;T|X,Y|/60; 

7  til  NKXT  X 

71  S  PRINT 
770  NKXT  Y 

77S  PRINT  TAR  (»)  !*-- - - " ; 

7  10  PRINT  TAIl(SH)  ;" - r - * 

7  IS  PRINT 

740  PRINT  “ARK  TIIKRF  ANY  CHANGES"; 

74  S  INPUT  C$ 

/SO  IF  CJ  I  1  ,  l|  ="  Y"  Til  KN  1000 
7SS  RETURN 

1  000  i  RKM  MARK  CHANGES 

loos  print  "which  solvent 

10  to  INPUT  .1 

lots  PRINT  “TO-*  ;V|.I|  /*>0;  *  ’I’ll  K  NKW  VAI.III-:  I  S“  ; 

10  70  INPUT  V|.!| 

l o 2 s  i.rr  v|.i|=v|.i|  »r»o 

t  010  PRINT  “  W I  I'T'I  =";W|J,l|;"  THE  NKW  VAI.UK  IK"; 

10 JS  INPUT  A 
1040  FOR  x-l  TO  M2 
1  (14  S  I, FT  W| J ,X | *A 

1OS0  NIXT  X 

10SS  PRINT  “WHICH  PKAK  NEEDS  TO  OK  CHANCED  I1-";M2;")"; 

1  Of)  0  INPUT  I 

1  06S  PRINT  "PKAK  “  ;  I  ;  “II  A*i  A  K  KTKNTI  ON  TIME  OF  "  ;  T’|.l  ,  I  )  /6  0;  "  TIIK  NEW  VALUE  IS" 
1  070  INPUT  T | J  ,  I  | 

TO 7 S  I. FT  T|.)  ,  I  |  -T|J  ,  I  |  *60 
1  ORO  I.  FT  K|.l  ,  I  I  =  <T  1  J  ,  II /V|  J|  )  -1 
10RS  PRINT 

1 0 10  PRINT  "ANY  MORF.  CHANGES  FOR  SOLVENT  "  ;.l  ; 

10  RS  INPUT  CS 

1100  IF  C?  J 1 , 1 1 *" Y"  THFN  105S 

1  10S  PRINT  "ANY  CHANGES  FOR  ANOTHER  SOLVENT"; 

11 1U  INPUT  i> 

Ills  ii  i  <  1  ,  i  |  «"Y*  THEN  100S 

117  0  (to sun  6io 

1  t  2 S  G.7T0  12  S 

A 000  REM  RESET  POINTERS  •** 

sons  pfah  »i , l 

6010  PEAn  #7,1 

60  IS  RE  An  » 3,  l 

6  0 20  RE  An  #4,1 

6  0  2S  RE  AO  * S, 1 

6030  RE An  16,1 

6  ’3S  RE  An  #7,1 

6040  PFTUHN 

qqog  rrjD 
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10  COM  A (2 1 ,T| 7, 35) 

100  REM  A************************************************* 

105  REM  ********  *LIQU  10  CHROMATOGRAPHIC  OPT  IKIZATICN*  ***  *» 

110  REM  ***********  XLOB6B  * *R. RADKE  /  J. CARVER  ********** 

115  REM  *  »*  ********************  »»AUC;/SEI>'r  J981*  *********** 

120  REM  ************************************************** 

125  CMRS  27, ES 

130  DIM  E$|l) , Y $ (101 , ZS | 26 J ,B(35) ,C|35) 

135  LET  Y$«" 0123 1567 89” 

140  LET  Z$="  ,  s  AM" 

145  LET  Mtlj»-M13)»M|51»M|71=M181-M(101-M|121-31 
150  LET  M I  4 1 =M l 6  J -M | 9 1 -M  J 1 1 ] =  3 0 
155  LET  M (2 1=28 

160  FILES  FIRST: : 19, SECOND: :19, THIRD: : 19 
165  FILES  FOURTH: : 19, FIFTH :: 19, SIXTH: : 19 
170  FILES  SEVNT1I:  :19,RESOLV:  :19, HOLDER:  :19 
175  READ  #8,1 

180  FOR  11=0  TO  100  STEP  2 

185  FOR  J1=0  TO  100-11  STEP  2 

190  PRINT  18; 1. 000008*38 

195  NEXT  J1 

200  NEXT  11 

2000  REM  ***  ORM  SUBROUTINE  *** 

2  00  5  OUTDVC  "T*"  ,E 
2010  LET  Ml=A  1 1 j 
2015  LET  M2=  A ( 2 | 

2020  GOSUB  6000 

2025  FOR  1=1  TO  M2-1 

2030  PRINT  “WORK ING I “ 

2035  FOR  K= I +1  TO  M2 

2040  GOSUB  6100 

2045  LET  V=0 

2050  GOSUB  3000 

2055  NEXT  K 

2060  NEXT  I 
2065  OUTDVC  "L2",E 
2070  PRINT  E$"U" 

2075  PRINT 
2080  PRINT 
2035  GOSUB  5000 

2090  PRINT  TAB (20) "HPLC  DATA  ANALYSIS" 

2095  PRINT  TAB (20) "RUN  "  ;  S  l  2 1 ; D$ [S { 3 ) - 3 ,S ( 3 ) ) ;  S 1 4 1 

2100  GOSUB  3300 

2105  GOSUB  3500 

2110  OUTDVC  "T*",E 

2115  PRINT 

2120  PRINT  "DO  YOU  WANT  TO  PLOT  LINES  OF  RESOLUTION?  (Y/N)"; 

2125  INPUT  C$ 

2130  IF  C$  f 1 , 1 ) I" Y "  THEN  STOP 
2135  PRINT 

2140  PRINT  "WHAT  IS  THE  MINIMUM  RESOLUTION  DESIRED  "; 

2145  INPUT  R9 

2150  PRINT  "PEAK  PAIRS  ALWAYS  RESOLVED  GREATER  THAN  ";R9;"WILL  BE  IGNORED 
2155  OUTOVC  " Cl " , E 
2160  PRINT  E$"U" 

2165  LET  Z  1=2 

2 ’70  REM  ***  GET  Y  AND  BETA  VALUES  *** 

2175  GOSUB  6000 

2180  FOR  1=1  TO  M2- 1 

2185  FOR  K= I *1  TO  M2 
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2190  GOSUH  6100 
2195  PRINT 
2200  PRINT 
2205  PRINT 

2210  REM  **♦  CHECK  FOR  R9.-R9,  AND  0  »** 
2215  OUTDVC  "T*",E 
2220  PRINT  E$"H",E$"J" 

2225  PRINT  "  FOR  PEAK  PAIR  ";Ij“-  "}K 
2230  LET  R7«Y«0 

2235  PRINT  "CHECKING  FOR  RESOLUTION  OF  ";Y 

2240  GOSUB  5500 

2245  IF  R7*l  THEN  2295 

2250  LET  Y*R9 

2255  PRINT  "CHECKING  FOR  RESOLUTION  OF  " s Y 

2260  GOSUB  5500 

2265  IF  R7»l  THEN  2295 

2270  LET  Y»-R9 

2275  PRINT  "CHECKING  FOR  RESOLUTION  OF  ";Y 

2280  GOSUB  5500 

2285  IF  R7«l  THEN  2295 

2290  IF  Zl-3  THEN  GOSUB  5700 

2295  RFM  **»  PRINT  HEADING  *** 

2300  PRINT  E$"b* 

2305  OUTDVC  "Cl" , E 

2310  IF  R7*l  THEN  PRINT  ES"U" 

2315  PRINT 
2320  PRINT 

2325  PRINT  "  FOR  PEAK  PAIR  "jlj"-  "}K 
2330  PRINT 

2335  PRINT  "  BETA  COEFFICIENTS  (1-7)“ 

2340  PRINT  Bl;B2»B3;B4;85j B6»B7 
2345  PRINT 
2350  PRINT 

2355  REM  ***  CHECK  POINTS  *‘* 

2360  PRINT  "  MEOH " , " ACN" , "THF" * " Y" 

2365  PRINT  " - 

2370  LET  X1-X2-.16 
2375  LET  X3*.68 
2380  GOSUB  5800 
2385  LET  X1-X3-.16 
2390  LET  X2-.68 
2395  GOSUB  5800 
2400  LET  XI  *  . 68 
2405  LET  X2-X3-.16 
2410  GOSUB  5800 
2415  LET  X1-X2-0 
2420  LET  X3*l 
2425  GOSUB  5800 
2430  LET  X1-X3-0 
2435  LET  X2 «1 
2440  GOSUB  5800 
2445  LET  Xl*l 
2450  LET  X2*X3»0 
2455  GOSUB  5800 

2460  PRINT  " - - - 

2465  PRINT 
2470  PRINT 

2475  IF  R7«l  THEN  LET  Zl-0 

2480  LEJT  Z1*ZU1 

2485  IP  R7«l  THEN  GOSUB  4000 
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2490 
2495 
2500 
2505 
2510 
2515 
2520 
2525 
2530 
2535 
2540 
3000 
300  5 
3010 
3015 
3020 
3025 
3030 
3035 
3040 
3045 
3050 
3055 
3060 
3065 
3070 
3075 
3080 
3085 
3090 
3095 
3100 
3105 
3110 
3115 
3120 
3125 
3130 
3300 
3305 
3310 
3315 
3320 
3325 
3330 
3335 
3340 
334  5 
3350 
3355 
3  360 
3  365 
3370 
3375 
3500 
3505 
3510 
3515 
3520 
3525 


IF  R7»l  THEN  PRINT  E$"& p5D" E$"& p5uOC" 

IF  R7I1  THEN  PRINT  "PEAK  PAIR  ";I;“-  ";K;"IS  ALWAYS  RESOLVED  "  ;  R9 
IF  R7=l  THEN  GOSUB  5700 
NEXT  K 
NEXT  I 

OUTDVC  *T*",E  . 

PRINT  E$"b" 

PRINT  "DO  YOU  WANT  TO  PLOT  A  DIFFERENT  RESOLUTION  (Y/N)"; 

INPUT  CS 

IF  C$(1,1)="Y"  THEN  2140 
STOP 

REM  •**  CALCULATES  BEST  SEPARATION  SOLVENT  *** 

READ  1 8 , 1 
READ  *9,1 

FOR  11=0  TO  100  STEP  2 

FOR  J1=0  TO  100-11  STEP  2 
LET  Xl-Il/100 
LET  X2=J 1/100 
LET  X3=1-X1-X2 

LET  Y=B1*X1+B2*X2+B3*X3+B4*X1*X2+B5*X1*X3+B6*X2*X3+B7*X1*X2*X3 
READ  18; Y1 

IF  AHS(Y)  <=  Y1  THEN  LET  Y1=ABS(Y) 

PRINT  #9; Y1 
IF  Y1  <«  V  THEN  3080 
LET  V*Yl 
LET  X6=I 1 
LET  X7  =J  1 
NEXT  J1 
NEXT  II 
READ  *8,1 
READ  *9,1 

FOR  11=0  TO  100  STEP  2 

FOR  J 1=0  TO  100-11  STEP  2 
READ  19 ; Y 
PRINT  1 8 ; Y 
NEXT  J1 
NEXT  II 
RETURN 

REM  •**  PRINT  BEST  SEPARATION  *** 

LET  XI-X6/100 

LET  X2=  X7/1 00 

LET  X3-1-X1-X2 

PRINT 

PRINT 

PRINT 

PRINT  TAB ( 15) ; "THE  BEST  COMPOSITION  IS:" 

PRINT 

PRINT  TAB (10) ; "MEOH "  ;TAB ( 2  5)  ;"ACN“;TAB (40) ; "THF" 

PRINT  TAB (1 1) ;X6;"%";TAB( 25) ;X7;"%";TAB (41) ; (100-X6-X7) 

PRINT 

PRINT  TAB (5) ; "THE  WORST  RESOLUTION  AT  THIS  POINT  IS  "jV 

PRINT 

PRINT 

RETURN 

REM  ***  PREDICT  ORDER  OF  ELUTION  *,** 

GOSUB  6000 
LET  I*C(1 ) =1 
LET  W=  20 
LET  B ( 1 1 *T 1 1 , 1 ) 

FOR  K=  2  TO  M2 
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3530  LET  C[KJ=K 
3535  GOSUB  6100 

3540  LET  Y»B1*X1+B2*X2+B3*X3+B4*X1*X2+B5*X1*X3+B6*X2*X3+B7*X1*X2*X3 
3545  LET  B  lK|  «W*Y+B  1 1) 

3550  NEXT  X 

3555  REM  ***  SORT  *** 

3560  FOR  1=1  TO  M2-1 

3565  FOR  J»I+1  TO  M2 

3570  IF  B[J|  >=  B ( I )  THEN  3605 

3575  LET  B1=B[I] 

3580  LET  Cl =C ( I  1 

3585  LET  C|I1*C[J] 

3590  LET  CiJi-Cl 

3595  LET  B  ( I ] “B [J] 

3600  LET  BJJ]»B1 

3605  NEXT  J 

3610  NEXT  I 

3615  PRINT  TAB(15) ;"PREDICTED  ORDER  OF  ELUTION* 

3620  PRINT 

3625  PRINT  TAB{10) ; "PEAK  (  *, "RELATIVE  RETENTION  TIME" 

3630  PRINT 

3635  FOR  1=1  TO  M2 

3640  PRINT  TAB (10) ;C( 1 1 ;TAB( 35) ;B [I )/6  0 

3645  NEXT  I 

3650  RETURN 

4000  REM  PLOT  GRID 

4001  OUTDVC  "C1",E 

4002  PRINT 

4003  PRINT 

4004  PRINT 

4005  PRINT 

4007  OUTDVC  "T*",E 

4010  DIM  I  (21  )  ,  J  ( 21  ]  ,L[ 21]  ,M(21],N[21),0[21) 

4030  LET  G1 *10 

4040  LET  G2M00 

4050  LET  Kl=ll 

4060  LET  D=*G2/  (Kl-1 ) 

4070  FOR  11=1  TO  K 1 
4080  LET  A=»  1 1-1 

4090  LET  I|I1J=INT{G1+(0*G2+(1)*A*D)/1) 

4100  LET  J[Il]-INT(Gl+(2*G2+(-l)*A*D)/2) 

4110  LET  L(Il]=INT(Gl+(l*G2+(-l)*A*D)/2) 

4120  LET  M[Il]=INT(Gl+(0*G2+(0) * A*D \ *SQR (3 ) /2 ) 

4130  LET  N(Ill=INT(Gl+(0*G2+(l)*A*D)*SQR(3)/2) 

4140  LET  0[I1)=INT(G1+ (l*d2+(-l) *A*D)*SQR(3)/2) 

4150  NEXT  II 

4160  REM  PRINT  GRID 

4170  PRINT  E$"H"E$"J"E$"*tn2a  0,0  jibZ"E$"*dachlntZ" 

4180  LET  M=  K1  -1 
4190  FOR  N= 1  TO  M 
4200  LET  L=N  +  1 
4210  LET  J»K1-N 

4220  PRINT  E$"*pah"I |L] ,M[L)" Z"E$"*pbh"L(J] ,0|J)"Z" 

4230  PRINT  E$"*pah"J[L] ,N[L]"Z"E$"*pbh"I(J) ,M(J]"Z" 

4240  PRINT  E$"*pah"L(L] ,0 (L ] " Z" E$ "*pbh " J [ J J ,N(J)"Z“ 

4250  NEXT  N 

4260  PRINT  E$ "*d 0 , OokSTHF" 

4270  PRINT  ES"*d422, OOACN " 

4280  PRINT  E$"*d230, 350OMEOH" 

4290  PRINT  E$"*dlT" 
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4  300  REM  ,  ***  PLOT  LINES  OF  RESOLUTION  *** 

4302  DIM  XflOll ,Y[101] ,PII01] ,Q[101) 

4303  LET  Y “P.9 

4304  GOSUB  4310 

4305  LET  Y  —  R9 

4306  GOSUB  4310 

4307  LET  Y»0 

4308  GOSUB  4310 

4309  RETURN 

4310  LET  N*M“1 

4311  FOR  X=0  TO  100 

4312  LET  X1=X/100 

4314  LET  Cl=B3+(Bl-B3+B5-b5*Xl)*Xl 

4316  LET  C2=B2-B3+B6+(B4-B5-B6+B7-B7*Xl)*Xl 

4318  LET  C3=0-B6-B7*Xl 

4320  LET  A  *C3 

4322  LET  B=C2 

4324  LET  C=C1-Y 

4326  LET  D=»B*B-4*A*C 

4328  IF  D<0  THEN  4370 

4330  LET  S=SQR(D) 

4332  LET  R1 = (0-B +S) / ( 2*A) 

4334  LET  R2 = (0-B-S ) / ( 2* A) 

4336  IF  R1<0  THEN  4352 

4338  IF  Rl> 1  THEN  4352 

4340  LET  01  =X1  * R1 

4342  IF  Q1>1  THEN  4352 

4344  LET  X2-1-Q1 

4346  LET  X3=R1 

4348  LET  X( N ) =Gl* (X3+X 1/2 ) *G2 

4350  LET  Y[N]=G1+-X1*G2*  SQR  (  3 )  /  2 

4351  LET  N  =N+ 1 

4352  IF  P2C0  THEN  4370 

4354  IF  R2>1  THEN  4370 

4356  LET  Q2=X1+R2 

4358  IF  Q2>1  THEN  4370 

4360  LET  X2=l-Q2 

4362  LET  X3  =  R2 

4364  LET  P (M ) »Gl+ 'X 3+X 1/2 ) *G2 

4366  LET  Q [M] =G1+X1*G2*SQR ( 3) /2 

4368  LET  M-M+l 

4370  NEXT  X 

4372  IF  N< 3  THEN  4392 

4374  LET  Y8  =  Y [ 1 | 

4376  LET  X8  =  X| 1) 

4378  FOR  19=2  TO  N-l 

4380  LET  X9=X [ 1 9 1 

4382  LET  Y9  =Y ( 19  I 

4384  PRINT  E$"*pah"X8, Y8"Z"E$"*pbh"X9, Y9"Z" 

4386  LET  X8=X 9 

4388  LET  vs=y9 

4390  NEXT  19 


4391 

PPINT 

E$"*dsZ";Y 

4392 

PRINT 

F,S"*(ltZ" 

4393 

IF  M<3 

THEN  RETURN 

4394 

•  LET  Y8 

=  0(1) 

4396 

LET  X8 

=  P  ( 1  1 

4398 

FOR 

19=2  TO  M-l 

4400 

LET 

X9=P ( I 9 ! 

4402 

LET 

Y9=Q ( I 9 ) 

.  ! 
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4404  PRINT  F$"»pah‘X8,y8"2"ES"*pbh"X9,V9"Z" 

4406  LET  XP-X9 
4408  LET  YP-«Y9 
4410  NEXT  19 

4412  PRINT  E$"*dsZ" 

4413  PRINT  Y 

4414  PRINT  E$"*dtZ" 

4416  RETURN 

5000  REM  ***  GET  DATE  AND  TIME  *** 

5010  DIM  D$[48] 

5020  LET  D$*"JAN  FEB  MAR  APR  MAY  JUNEJULYAUG  SEPTOCT  NOV  DEC  " 

5021  LET  Z2“TIM (3) 

5022  LET  Z3«9 

5023  LET  Z4-12 

5024  IF  Z2/400-INT(Z2/400)  LET  M 1 2 1 -29 

5025  IF  Z2/100IINT(Z2/100)  AND  Z2/4-INT ( Z2/t >  LET  M[2)«29 

5026  LET  Z  2  =T  I M  (  2 ) 

5027  FOR  Z=1  TO  12 

5028  LET  Z2=Z2-M{ Z] 

5029  IF  Z2  <«  0  THEN  5031 

5030  NEXT  Z 

5031  LET  S [2 ) »«2+M ( Z ] 

5040  LET  S [ 4 ) »TI M ( 3) 

5060  LET  S[3]«Z*4 
5070  RETURN 

5500  REM  ***  CALCULATES  PREDICTED  RESOLUTION  VALUES  *** 

5502  PRINT  E$“c" 

5504  FOR  X“0  TO  100  STEP  2 

5506  LET  X1*X/100 

5508  LET  C1»B3+(B1-B3+B5-35*X1)*X1 

5510  LET  C2»B2-B3+B6  +  (B4-B5-B6  +  B7-B7* Xl) * XI 

5512  LET  C3»0-B6-B7*Xl 

5514  LET  A=C3 

5516  LET  B=*C2 

5518  LET  C«C1-Y 

5520  LET  n-B“B-4*A*C 

5522  IF  0<0  THEN  5548 

5524  LET  S»SOR(D) 

5526  LET  R 1 = ( 0-B+S) / ( 2* A) 

5528  LET  R2 • ( 0-B-S) / ( 2* A) 

5530  IF  R1<0  OR  R1 > 1  THEN  5538 
5532  LET  01=X1+R1 

5534  IF  QI>1  THEN  5538 

5536  LET  F7=l 

5538  IF  R2 < 0  OR  R2>1  THEN  5546 

5540  LET  Q2*X1+R2 

5542  IF  Q2>1  THEN  5546 

5544  LET  R7=l 

5546  IF  R7=l  THEN  RETURN 

5548  NEXT  X 

5550  RETURN 

5700  REM  ***  MEW  PAGE  *** 

5702  OUTDVO  “Cl" , B 
5704  PRINT  FS’U” 

5706  LET  Zl=l 
5708  RETURN 

5K00  REM  ***  CALCU LATE  PREDICTED  RESOLUTION  VALUES  *** 

5802  LET  Y=nl*XU-B2*X2+B3*X3+U4»Xl»X2+n5*Xl*X3+B6*X2*X3+B7*Xl*X2*X3 
5804  PRINT  Xl*100,X2*100,X3*100,Y 
5806  RETURN 
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This  is  basically  a  translated  version  of  XLDB6A  and  B  but  with  many 
enhancements.  It  still  has  a  section  to  plot  lines  of  resolution  on  a 
Hewlett-Packard  graphics  terminal  with  a  graphics  printer  capable  of  a  RASTER 
Dump  but  it  is  not  limited  to  a  graphics  terminal.  Extensive  use  of  the 
Hewlett-Packard  graphics  package  "DGL"  has  been  made.  This  allows  the  user  to 
use  a  4  pen  graphics  plotter  for  output  from  a  nongraphics  terminal.  Figure  1 
of  report  text  is  typical  of  the  output.  This  greatly  simplifies  selection  of 
a  solvent.  Another  enhancement  is  that  the  output  gives  expected  retention 
time  of  named  peaks  for  any  selected  composition  and  the  minimum  resolution 
predicted  at  that  composition.  Table  3  of  report  text  is  an  example  of  this 
output. 
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CHARACTER  REPRESENTATION  OF  COLUMN  NAME 
CHARACTER  REPRESENTATION  OF  SCALE 


INITIALIZATION 

PROGRAM  LCMAIN 

COMMON  /ONE/NUMP  EAK ,DATE , BESTX1 , BESTX2  ,BESTX  3 .PEAKNAME 
COMMON  /TWO /BEST RES ,M INRES , WIDTH , KPRIME,PEAKRT , VOIDRT 
COMMON  /THREE/RESOLUTION 

COMMON  /FOUR/PEN1 X , PEN1 Y ,  PEN2  X , PEN2  Y,PEN3X,PEN3Y,  PEN4X  ,PEN4 Y 
COMMON  /F I VE/PEN 1PT, PEN 2PT , PEN3 PT, P EN4PT , I  SAMPLE, I COLUMN, 

1  I  MEOH, I  THF , I _ACN , I SC  ALE 
COMMON  /StX/XARRAY  ,Y ARRAY ,  PARRA  Y,  QARRA If 
COMMON  ESC , FF , LU 

INTEGER  NUMPEAK,LU,DATE(15)  , ORDER (35) ,CONTRL, IERR , PENlPT , 

1  PEM2PT ,PEN  3PT , PEN4PT , ISAMPLE (8) ,1 COLUMN ( 8 ) , I  MEOH(2),I  THF (2) , 

1  I  ACN ( 2 ) , I SCALE ( 3) ,PAGE , PLOT, XARRAY (101 ) ,Y ARRAY (101) , 

1  PARR AY ( 101 ) ,Q ARRAY (101) 

REAL  COF ( 7 )  ,VO IORT (7)  ,WIDTH(7,35) , KPRIME (7, 35) ,PEAKRT (7 , 35) , 

1  RESWANT, RESOLUTION (7, 35, 35)  .MINRES (0 : 50 ,0s  50)  ,BESTRES .TIME ( 35) , 
1  B  E  ST  X 1 , Q  E  ST  X2,BESTX3,XB(0:10),YB(0:10) ,XL (0 : 10) ,YL(0s 10) , 

1  XR (0 : 10 ) , Y  R ( 0 : 10 )  ,PEN1X ( 1 300) , PENl Y ( 1300)  ,PEN2X(1000)  , 

1  PF.N2Y  ( 1000  )  ,PEN3X(  1000)  ,PFN3Y  (1000  ),  PEN4X  ( 1000 )  ,  PEN4Y  ( 1000 )  , 

1  WIDE, HEIGHT ,VMAX,VMIN 

CHARACTER  PEAKNAME ( 35) *  15 , SOLV { 7 ) * 9 , F I LEI* 14 , C  MEOH*4,C  ACN»4, 

1  C  THF*4,CSAMPLE« 16,CCOLUMN*16,CSCALE»6 
CHARACTER  *1  ESC, FF, QUEST 
DIMENSION  LUBUF (275) 

EQUIVALENCE  (C  MEOH, I  MEOH),(C  THF , I_THF ) , (C  ACN  ,  I  ACN), 

1  (CSAMPLE, ISAMPLE) , (CCOLUMN , I COLUMN ) , ( CSCALETl SCA LE) 

DATA  SOLV  (  1 )  ,  SOLV  (2)  ,SOLV  (3)  /  'MEOH'  ,  'ACtl'  ,  'TIIF  '/ 

DATA  SOLV  (  4 )  .SOLV  (5)/  'MEOIUACN  ',  'MEOIUTHF  '/, 

DATA  SOLV ( 6) ,SOLV (7 )/  ' ACN6THF  ' , 'ALL  THREE'/ 

ESC=CHAR ( 27 ) 

FF=CIIAR  (12) 

PEN 1 PT=0 
PF.N2PT=0 
PEN  3PT=0 
PCN4PT=0 

CHANGE  THE  FILE  TO  FIT  YOUR  SYSTEM  ! 

FILE1  (7: ) ** ; :JC: 3: 9'  1  PREPARE  FILE  TO  BE  ON  LU  JC, TYPE  3 

CALL  LGHUF  ( LUiJ'JF  ,275) 

WRITE  (1,290)  ESC, ESC  !  CLEAR  SCREEN  AND  PRINT  QUESTION 


PEAKNAME (35) *15 

SOLV (7 )*9 

QUEST* 1 

FF*  1 

E5C*1 

FILE1*14 

C  MEOH* 4 

c"acn*< 

C— THF*  4 
C3AMPLE* 
CCOLUMN* 16 
CSCAL  E*  6 


C-4 


123  READ  (1,294)  1 COLUMN 

124  WRITE (1, * ( “  ENTER  THE  SAMPLE  NAME:  “) ') 

125  READ  (1,294)  ISAMPLE 

126  WRITE ( 1 ,  '  (  “  DO  YOU  WANT  TO  USE  AN  OLD  DATA  PILE?  ") *: 

127  READ  (1,391)  QUEST 

128  IP  (QUEST. NE.  'Y')  THEN 

129  WRITE  (1,291) 

130  READ  (1,293)  PILE1(:6) 

131  OPEN  <333,FILE-PILE1,IOSTAT-IOS1,ERR-910,STATUS-'HEW') 

132  ELSE 

133  WRITE  (1,292) 

134  READ  (1,293)  PILE1(:6) 

135  OPEN  (333,PILE>FILE1, IOSTAT-IOS1 , ERR- 910 .STATUS*  'OLD ' ) 

136  REWIND  (333) 

137  READ  (333)  PEAKNAME 

138  DO  105  1-1,7 

139  READ  (333)  VO IDRT ( I ) , (PEAKRT (I , J) .WIDTH ( I , J ) , J-l , 35) 

140  105  CONTINUE 

141  URITE  (1,289) 

142  READ  * , MUMP  EAR 

143  IF  (NUfPEAK.EQ.O)  STOP 

144  GOTO  225 

145  ENDIF 

146  * 

147  **»*  GET  PEAK  DATA 

148  * 

149  200  WRITE  (1,289) 


150  READ  *,NUPPEAK 

151  IF  (NUMPEAK . EQ. 0)  STOP 

152  IF  (NUMPEAK.LT. 2. OR. NUMPEAK .GT. 35)  GOTO  200 

153  DO  220  1-1,7 

154  WRITE  (1,280)  SOLV(I) 

155  READ  * .VOIDRT ( I ) 

156  IF  (VOIDRT (I ) .EQ. 0)  STOP 

157  IF  (I.LT.2)  THEN 

158  WRITE  (1,281) 

159  DO  205  J-l, NUMPEAK 

160  WRITE  (1,282)  J 

161  READ  * , PEA KRT ( I , J ) ,W I DTH ( I , J ) , PEAKNAME (J ) 

162  205  CONTINUE 

163  ELSE 

164  210  WRITE  (1,283) 

165  DO  215  J-), NUMPEAK 

166  WRITE  (  1,284)  PEA  KNAME ( J) 

167  READ  *, PEAKRT ( I ,J) ,WIDTH(I,J) 

168  215  CONTINUE 

169  ENDIF 

170  220  CONTINUE 

171  * 

172  *****  PRINT  INPUT 

173  * 

174  225  WRITE  (1,285)  VOIDRT 

175  DO  230  I -1, NUMPEAK 

176  WRITE  (1,286)  I , PEAKNAME ( I ) , ( PEAKRT ( J , I ) , J-l , 7) 

177  WRITE  (1,287)  (WIDTH  (J, I), J-l, 7) 

178  230  CONTINUE 

179  WRITE  (1,288) 

180  280  FORMAT  (/.'ENTER  THE  RETENTION  TIME  (MIN)  FOR  SOLVENT  *,A9, 

181  1'  <0-EXIT> ' ) 

182  201  FORMAT  (/.'ENTER  THE  RETENTION  TIME  (MIN),  PEAK  WIDTH  (SEC), 

183  1'  PEAK  NAME',/, 'PEAK  t  RT  WIDTH  NAME') 


C-5 


•>  i  "V 


184  282 

185  283 

186 

18?  284 

188  285 

189 

190 

191 

192 

193  286 

194  287 

195 

196  288 

197 

198  289 

199  290 

200 

201  291 

202  292 

203  293 

204  294 

205  295 

206  • 

207  *•••• 

208  * 

209  300 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220  305 

221 

222  310 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
2  38 

239 

240 

241 

242 

243  315 


FORMAT  (T2,I2,3X,“  *) 

FORMAT  (/, 'ENTER  RETENTION  TIME  (MIN),  AND  PEAR  WIDTH  (SEC)',/, 

1'  PEAK  RT  WIDTH') 

PORMAT  (T2 ,A1 5, IX  ,*  *) 

FORMAT  (/,T48,  'ACN6  *7t57  , 'TIIFfc  ' ,T66,  'ACN6 *,T74 ,  'MEOfU ' ,/»T21 , 
l'MEOtl  ACN  TllF  MEOII  MEOII  T1IF  ACIUTHF ',/, 

1T2  2, '1 ',T31, '2  ',T40, '3',T49, *4',T58,  '5',T67, '6 ' ,T76 , '7 ' ,/ ,T19, 

1  ' - - - — - — - — - - - -  */ 1 

1T5, 'SOLVENT  FRONT ' ,T21 ,F5 . 2 , 6F9 .2 ,/) 

FORMAT  (I3,X1,A15,T21,F5.2,6F9.2>  ,  .  , 

FORMAT  (T22, '( ',13,  ') ',T31,  '( ',13,  ') ',T40,'( ',13,  )  ,T49,,  (  ,13, 
1')  ',T58,  '(',13,  ')  '  ,T67,  '(',13,  ')  \T76,'(',13,  ')  ') 

FOfflAT  (T19,  ' -  » 

1' . ') 

FORMAT  ('ENTER  THE  NUMBER  OF  PEAKS  PER  SOLVCUT.  (2-35)  <0«EXIT>  ) 
FORMAT  (Al, 'II', Al, 'J ',//,T30,'HPLC  OPTIMIZATION',//, 

1'  ENTER  THE  COLUMN  NAME:  ') 

FORMAT  ('  ENTER  THE  NAME  OF  THE  NEW  FILE.  ') 

FORMAT  ('  ENTER  THE  NAME  OF  THE  OLD  FILE.  ”') 

FORMAT  (A6) 

FORMAT  (8A2) 

FORMAT  (Al, 'll'.Al, 'J') 


CHANGE  SECTION 


WRITE  U,'(*  ARE  THERE  ANY  CHANGES?  *)') 

READ  (1,391)  QUEST 
IF  (QUEST. NE.  'Y ' )  GOTO  320 

WRITE  (1,'("  DO  YOU  WANT  TO  CHANGE  A  NAME?  _*) ') 

READ  (1,391)  QUEST 

IF  (QUEST. NE. 'Y ')  GOTO  305 

WRITE  (1,'C  ENTER  THE  PEAK  I  *)') 

READ  «,I 

WRITE  (1,394)  I  , PEAKNAME(I ) 

READ  * , PEAKNAME ( I ) 

GOTO  315 

WRITE  (1,'C  WHICH  SOLVENT?  _*)') 

ABAC)  *  I 

WRITE  ( 1 , ' ( "WHAT  DO  YOU  WANT  CHANGED?  ( VOIDRT  ,PEAKRT , OR  WIDTH)*)') 
READ  (1,391)  QUEST 
IF  (QUE'i EQ.  'V'  )  THEN 

WRIT'  (1,385)  VOIDRT ( I ) 

READ  * , VOI DRT ( I ) 

ELSEIF  (QUEST. EQ.  'P ' )  THEN 

WRITE  (1,'(*  WHICH  PEAK  I  ?  *)') 

READ  *,J 

WRITE  (1,386)  PEAKNAME ( J) , PEA  KRT (I , J ) 

READ  *,PEAKRT(I,J) 

ELSEIF  (QUEST  .EQ. 'w')  THEN 

WRITE  (1,'C  WHICH  PEAK  I  ?  *)') 

READ  *,J 

WRITE  (  1,387)  PEAKNAME  (J  )  ,WIDTII ( I ,  J ) 

READ  *, WIDTH  (I,J) 

ELSE 
GOTO310 
END  IF 

WRITE  (1,'(*  ANY  OTHER  CHANCES  FOR  THIS  SOLVENT?  ")') 

READ  (1,391)  QUEST 

IF  (QUEST. EQ. 'Y ')  GOTO  310 

WRITE  ( 1 ,  '  ( "ANY  OTHER  CHANGES?  ")') 

READ  (1,391)  QUEST 


06 


245 

246 

247 

248 

320 

249 

250 

251 

2  52 

325 

253 

254 

255 

2S6 

257 

258 

2S9 

260 

261 

330 

262 

263 

264 

265 

266 

267 

268 

335 

269 

270 

271 

272 

273 

274 

275 

276 

340 

277 

345 

278 

279 

280 

281 

282 

283 

350 

284 

355 

285 

360 

286 

287 

288 

28  9 

290 

365 

291 

370 

292 

375 

29  3 

294 

C 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

380 

IP  (QUEST. EQ. 'Y ')  GOTO  305 
GOTO  225 

WRITE  T1IE  DATA  TO  TOE  DISC  PILE 
REWIND  (333) 

WRITE  (333)  PEAKNAME 
DO  325  I-  1,7 

WRITE  (333)  VOIDRT (I ) , (PEAKRT (I , J ) .WIDTH ( I , J ) ,J«1,35) 
CONTINUE 
CLOSE  (333) 

WRITE  <1,'( "DO  YOU  WANT  TO  ANALYSE  THE  ENTIRE  LIST?  <Y/N>  ") 

READ  (1,391)  QUESr 
IF  (QUEST. EQ. *Y')  GOTO  360 

WRITE  ( 1 , * ( " HOW  MANY  PEAKS  DO  YOU  WANT  TO  ANALYSE?  ")') 

READ  * ,NUHPEAK 
DO  335  I -1, NUMPEAK 

WRITE  (1, '("ENTER  A  PEAK  NUMBER.  *)') 

READ  *, ORDER (I) 

WRITE  (1,395)  ORDER ( 1 ) , PEAKNAME (ORDER ( I ) ) 

READ  (1,391)  QUEST 
IF  (QUEST. NE. 'Y')  THEN 

WRITE  (1, '("RE-ENTER  THE  PEAK  #.  ")') 

GOTO  330 
ENDIF 
CONTINUE 

DO  345  I«1,NUMPEAK-1 
DO  340  J-I+l  .NUMPEAK 

IF  (ORDER(J) .LT.ORDER(I))  THEN 
ORD1-ORDER (I ) 

ORDER(I) -ORDER (J) 

ORDER ( J ) -ORDl 
ENDIF 
CONTINUE 
CONTINUE 

DO  355  I  "1  ,  NUMPEAK 

PE:..r::.'«ME(i)  -peakname (order  (i)) 

DO  350  J-1,7 

PEAKRT (J , I ) «P EAKRT (J , ORDER ( I ) ) 

WIDTH (J.I) -WIDTH (J, ORDER (I ) ) 

CONTINUE 
CONTINUE 
DO  370  1-1,7 

VRT-VOIDRT(I)*6n 
DO  365  J-l  .NUMPEAK 

RT -PEAKRT (I , J ) *60 
KPRIME ( I , J ) - (RT-VRT) /VRT 
CONTINUE 
CONTINUE 
WRITE  (1,388) 

READ  * , LU 

CHANGE  THE  LU  NUMBERS  TO  FIT  YOUR  SYSTEM 
IF  (LU.EQ.6)  WRITE  (LU,389)  l  FORM  FEED  LINE  PRINTER 
IF  (LU.EQ.7.0K.LU.EQ. S3) 'WRITE  (LU,*)  FF  1  FF  9871  OR  9876 
IF  (LU. NE .1 .AND.LU. NE .6 .AND.LU .NE  .7.AND.LU .NE.53)  GOTO  375 
CALL  FT 1 ME (DATE) 

WRITE  (LU,39Q)  I  SKIP  3  LINES 

WRITE  (LU , 392)  DATE , CCOLUMN , CS AMPLE 
WRITE  (LU , 28  5)  VOIDRT 
DO  380  1-1, NUMPEAK 

WRITE  (LU, 286 )  I , PEAKN AME { I ) , (PEAKRT( J, I ) , J- 1 , 7) 

WRITE  ( LU , 2  8  7 )  (WI DTH ( J , I )  , J-l, 7 ) 

CONTINUE 
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I 


I 
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/ 


/  i 


i 


z  V 


306  WRITE  (LU ,288) 

307  385  FORMAT  ("  VOIDRT-  ",F5.2,"  NEW  VALUE  IS  TO  BE  •) 

308  386  FORMAT  (T2,A15,“  RETENTION  TIME  »  ",F7.2,"  NEW  VAI UE  ",*  ") 

309  387  FORMAT  (T2,A15,"  PEAK  WIDTH  »  *,F7.2,"  NEW  VALUE  “") 

310  C  CHANGE  THE  LU  NUMBERS  TO  FIT  TOUR  SYSTEM 

311  388  FORMAT  (/,*  WHERE  DO  YOU  WANT  THE  OUTPUT  ?",/,T10,"l  -  THIS  TERM" 


s'  < 

312 

1" INAL  * ,/ ,T10 , *6  -  THE  LINE  PRINTER" ,/,T10 , "7  -  THE 

DAISY 

WHEEL  ", 

\ 

313 

1" PRINTER" ,/,T10," 53-  THE  THERMAL  PRINTER") 

% 

314 

389 

FORMAT  ("1") 

315 

390 

FORMAT  (/,/) 

s* 

316 

391 

FORMAT  (Al) 

» 

• 

| 

317 

392 

FORMAT  (T30  ,  'llPLC  OPTIMI Z Al ION ' ,//,T2 3, 15A2 ,//,T34 , 

'INPUT 

DATA',/, 

318 

1T20,  'COLUMN  s ' ,A16 ,T45 ,  'SAMPLE  s ' ,A16,//,T25, ' (DATA 

IS  CORRECTED  , 

V.* 

319 

1'  TO  1 .00  ml/min) ') 

320  393  FORMAT  (A15) 

FORMAT  {'  PEAK  ',12,'  'Z  CURRENTLY  ',A15,/,  '  ENTER  THE  HEW 
1  'NAME  ') 

FORMAT  (7. 'PEAK  •  ',12,'  IS  CURRENTLY  ',A15,'  OK?  <Y/N>  ') 


321  394 

322 

323  395 

324  • 

325  ***** 

326  * 

327 

328 

329 

330 
3  31 

332  400 

333  405 

334  410 

335 

336 

337 

338  415 

339 

340 

341 

342 

343 
3*4 


CALCULATE  RESOLUTION 


• 

345 

346 

» 

■ 

347 

1 

348 

420 

r 

349 

425 

350 

4  30 

v‘ 

351 

• 

j** 

352 

•  »* 

\: 

353 

* 

354 

i 

355 

C 

356 

h 

357 

L 

358 

r- 

359 

360 

[* 

361 

k 

362 

435 

J 

363 

P  , 

364 

r 

365 

365 

K 

h 

DO  410  1-1,7 

DO  405  J-1,HUMPEAK-1 
DO  400  K-J  +  l , NUMPEAK 

RESOLUTION  ( I , J ,  K)  -  (PEAKRT  (I ,  K) *60..-PEAKRT(  I , J )  *60 . )/ 

I  (  (WIDTH  ( I,  K)  +WIDTIl(I,J))*0.5) 

CONTINUE 

CONTINUE 

CONTINUE 

WRITE  (1,'(*  DO  YOU  WANT  TO  CALCULATE  THE  COF  ?  <Y/M>  ")') 

READ  (1,391)  QUEST 
IF  (QU..ST.EQ.  'Y ')  THEN 

WRITE  (1,'(*  ENTER  THE  DESIRED  RESOLUTION.  ")') 

READ  * , RESWANT 
IF  (RESWAHT.LT. 0)  GOTO  415 
DO  430  1-1,7 
COF ( I ) -0 

DO  425  J -1 ,HUMP EAK-1 
DO  420  K-J+l, NUMPEAK 

2 -RESOLUTION { I , J , K ) /RESWANT 
IF  (2.LE.0)  GOTO  420 
COF ( I ) -COF ( I ) +ALOG (Z ) 

CONTINUE 

CONTINUE 

CONTINUE 

OUTPUT 

DO  450  1-1,7 

CHANGE  THE  LU  NUMBERS  TO  FIT  YOUR  SYSTEM 

IF  (LU.EQ.6)  WRITE  (LU,389)  I  FORM  FEED  LINE  PRINTER 
IF  (LU.EQ.7.0R.LU.EQ.53)  WRITE  (LU,*)  FF  1  FF  9871  OR  9876 
WRITE  (L’J  ,  490)  SOLV  ( I )  ,  VOIDRT  (I ) 

DO  435  J-l, NUMPEAK 

WRITE  (LU , 491)  PEAKNAME (J ) , PEAKRT ( I , J ) , WIDTH ( I , J ) , 

1  KPRIME (I , J) 

CONTINUE 
WRITE  (LU, 492) 

DO  445  J-l ,NUMPEAK-1 
DO  440  K-J  +  l, NUMPEAK 

I  3- INDEX (PEAKNAME (J ) , '  ') 
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367 

368 

369 

370  440 

371  44S 

372 

373  450 

374 

375  490 

376 

377 

378  491 

379  492 

380  493 

381  494 

382 

383  * 

384 
185  • 

386  500 

387  505 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399  C 

400 

401 

402 

403  * 

404  *•*«< 

405  * 

406  910 

407  980 

408 


IF  (I3.GT. 15. OR. I3.LT. 1)  13-15 

WRITE  (LU  ,493)  PE.AKNAME  (J  )  ( :  1 3)  ,  PEAKNAME  ( K)  , 

1  RESOLUTION  (I ,  J , K) 

CONTINUE 

CONTINUE 

WRITE  (LU > 494)  SOLV (I) ,COF (I ) 

CONTINUE 

ENDIF 

FORMAT  (//,T5,'TIME  2ERO  rOR  SOLVENT  ',A9,'  IS*,F7.2, 

1'  MIN.',/,/,'  INGREDIENT' ,T20, 'RETENTION  TIME ' ,T37 , 'QAND  WIDTH* 
IT 50, 'CAPACITY',/) 

FORMAT  (T2,A15,T25,F5.2,T40,I3,T51,*6.2) 

FORMAT  .’//,T8,  '  FOR  PEAKS' ,T35, 'RESOLUTION ',/) 

FORMAT  (T3,Ri5,'-  ' ,A1 5 ,T36  ,F8 . 3 ) 

FORMAT  (//,'  CHROMATOGRAPHIC  OPTIMIZATION  FUNCTION  (COF)  FOR  ', 
1A9 , '  IS\F9.4) 


CALLS  TO  SUBROUTINES 

GO-0 

CALL  LCSORT  (GO) 

WRITE  (l, ‘("DO  YOU  WANT  A  PLOT  OF  THE  RESOLUTION  POINTS?  *)') 

READ  (1,391)  QUEST 

IF  (QUEST. EQ. 'Y')  CALL  LCPOINT 

WRITE  ( 1 , ' ("DO  YCU  WANT  A  PLOT  OF  THE  RESOLUTION  LINES?  “)') 
READ  (1,391)  QUEST 
IF  (QUEST. EQ.  'Y')  CALLLCLINE 

WRITE  ( 1 , ' ("DO  YOU  WANT  TO  CALC  FOR  A  DIFFERENT  SOLVENT?  ")') 

READ  (1,391)  QUEST 

GO^l 

IF  (QUEST. EQ. 'Y')  GOTO  505 
WRITE  (1,'(“  THE  END  I")') 

CHANGE  THE  LU  NUMBERS  TO  FIT  YOUR  SYSTEM 

IF  (LU.EQ.6)  WRITE  (LU.389)  I  FORMFEED  LINE  PRINTER 

IF  (LU .  EQ.  7  .  OR  .LU.  EQ.  53)  WRITE  (LU,*)  FF  1  FF  9071  OP.  5d76 

STOP 

******************  ERROR  SECTION  ****************************** 

WRITE  (1,780)  IOS1 

FOK1AT  (.  /, 'ERROR  IN  OPEN  STATEMENT.  IOSTAT  -  ',14) 

END 


* 


* 
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409  *  * 

410  A**********************  BEST  SEPARATION  SECTION  ********************** 

411  *  » 

412  $EMA  /TWO/,AHREE/,/FOUIV',/SIX/ 

413  SU3R0UTINE  LCSOKT  (GO) 

414  COMMON  /ONE/NUMPEAK  .DATE  ,  BESTX1 ,  BESTX2  ,BESTX3 ,  PEAKtJAME 

415  COMMON  /TWO /BESTRES .MINRES .WIDTH , KPRIME ,PEAKRT,VOIDRT 

4 IS  COMMON  AHRE  E/RESOLUTION 

41?  COMMON  /F I VE/PEN 1PT  t PEN 2PT.PEN3PT.PEN4PT, I SAMPLE. I COLUMN , 

418  1  1  MEOH,  I  TIIP, I  ACN ,  I  SC  ALE 

419  COMMON  ES?,FF,LU 

420  INTEGER  NUMPEAK ,LU, DATE  (15)  ,ORO  ER  ( 35)  ,CONTRL,  I  ERR,  PEN1PT, 

421  1  PEN2PT.PEN3PT ,  PEN4PT,  I  SAMPLE  (e)  ,  I  COLUMN  ( 3 )  ,  I  MEOH(2),I  TI1F  (2)  , 

422  1  I  ACN (2) , I SCALE ( 3) .PAGE , PLOT 

423  RERl  VOIDRT(7)  ,WIDTH(7,35),KPRIME(7,35)  .PEAKRT ( 7 , 35 ) , 

424  1  RESOLUTION ( 7,35,  35) , MIN RES (Ot  50, Ot  50) ,BESTRES,TIKE(35) , 

425  1  BESTXl, BESTX2 ,BE5TX3 

426  CHARACTER  PEAKNAME (  35)  *15 ,FILE1*14 ,C  MEOHM.C  ACN*4, 

427  1  C  THF*4 .CSAKPLE* 16 ,OCOLUMN*16,CSCALE*6 

428  CHARACTER  *1  ESC, FF, QUEST 

429  EQUIVALENCE  (C  MEOH, I  MEOH),(C  THF , I  THF),(C  ACN, I  ACN), 

430,  1  (CSAMPLE, ISAMPLE) , (CCOLUMN , I COLUMN ) 7 (CSCAL ETiSCALE) 


431 

IF 

(GO.GT.O)  GOTO  145 

432 

100 

DO 

110  1-0,50 

433 

DO  105  J-0,50-1 

434 

MINRES ( I , J ) - 1 . OE  38 

435 

105 

CONTINUE 

436 

110 

CONTINUE 

437 

WRITE  (1,388)  ESC, ESC  i  REVERSE  VIDIO  AND  PRINT  "WORKING” 

438 

DO 

130  1-1, NUMPEAK- 1 

439 

IF  (I.GT.l)  WRITE  (1,389)  ESC, I, ESC  !  REVERSE  VIDIO  AGAIN 

440 

DO  125  J-I+l, NUMPEAK 

441 

BETA1  -  RESOLUTION  (1,1, J) 

442 

BETA 2  -  RESOLUTION  (2,1, J) 

443 

BETA 3  -  RESOLUTION  (3,I,J) 

444 

BETA4  -  4.0 ‘RESOLUTION  (4,I,J)  -  2. 0* (BETA1+BETA2) 

445 

BETA 5  «  4. O'RESOLUTION  (5,I,J>  -  2 . 0* (BETA 1+BETA3) 

446 

BETA 6  -  4. O'RESOLUTION  (6,I,J)  -  2.0* (BETA2+BETA3) 

447 

BETA 7  -  27. O'RESOLUTION ( 7 , I , J )  -  12 . 0* (RESOLUTION ( 4, I ,J)  ♦ 

448 

1 

RESOLUTION (5,1 ,J) 'RESOLUTION (6, I , J) )+3. 0* (BETA1+BETA2+BETA3) 

449 

DC  120  11-0,100,2 

450 

DO  115  Jl-0, 100-11, 2 

451 

Xl-Il/100. 

452 

X2-J1/100. 

453 

X3-1. 0-X1-X2 

454 

Y-BETAI *X 1+BETA2*X2+BETA3*X3'BETA4*X1*X2+BETA5*X1*X3+ 

4  55 

1 

BETA6*X2*X3+BETA  7*X1*  X2*X3 

456 

Y-ABS (Y, 

457 

I 2-11/2 

458 

J2-J1/2 

459 

IF  (Y.LE.MINRES(I2,J2) )  MINRES ( I 2, J2) -Y 

460 

115 

CONTINUE 

461 

120 

CONTINUE 

462 

125 

CONTINUE 

463 

130 

CONTINUE 

464 

*** 

SELECT  THE  BEST  COMPOSITION  AND  ITS  WORST  RESOLUTION 

465 

WRITE  (1,380)  ESC, ESC  1  CLEAR  THE  SCREEN 

466 

BESTRES-0 

467 

DO 

140  1-0,50 

468 

DC  135  J-0 ,50-1 

469 

IF  (MINRES (I, J) .GT.BESTRES)  THEN 

C-10 


I  . 


v  w.1 


470  BESTRES-MINRES(I.J) 

471  BESTX1=I*2.0 

472  BESTX2»J*2. 0 

473  BESTX3-100.-BESTX1-BESTX2 

474  END IE 

475  135  CONTINUE 

476  140  CONTINUE 

477  * 

478  ***************  PRINT  BEST  SEPARATION  ********************** 

479  * 

480  C  CHANGE  THE  LU  NUMBERS  TO  FIT  YOUR  SYSTEM 

481  145  IE  (LU.EQ.6)  WRITE  (LU,381)  !  FORM  FEED  LINE  PRINTER 

482  IP  (LU.EQ.7.0R.LU.EQ. 53)  WRITE  (LU,*)  EE  !  FF  9871  OR  9876 

483  WRITE  (LU , 382 )  I  SKIP  3  LINES 

484  WRITE  (LU,385)  DATE ,CCOLUMN ,CSAMPLE, BESTX1 ,BESTX2  ,BESTX3 , BEST RES 

485  * 

486  ***************  PREDICT  ORDER  OF  ELUTION  ******************* 

487  * 

488  CI!0SX1»BESTX1 

489  CHOSX2=BESTX2 

490  CHOSX3-BESTX3 

491  RES-BESTRES 

492  WRITE  ( 1 , * ("WHAT  PLOW  RATE  ARE  YOU  INTERESTED  IN?  ")') 

493  READ  * , FLOW 

494  WRITE  ( i ,  * ("DO  YOU  WANT  TO  CHOSE  YOUR  OWN  COMPOSITION?  ")') 

495  READ  (1,383)  QUEST 

496  IF  (QUEST. EQ. 'Y')  THEN 

497  WRITE  (1,385)  DATE , CCO LIMN, CS AMPLE , BESTX1, DE5TX2 ,BESTX3 .BEST RES 

498  WRITE  (1, '("ENTER  COMPOSITION  YOU  WANT  TO  CALC  RES  TIMES  FOR")') 

499  READ  * .CHOSXl ,CHOSX2,C!iOSX3 

500  I»  INT (CHOSX1/2 . 0) 

501  J«INT (CHOSX2/2 . 0) 

502  RES«MINRES(I ,J) 

503  IF  (LU.EQ.6)  WHITE  (LU,301)  !  FORM  FEED  LINE  PRINTER 

504  IF  (LU.EQ.7.0R.LU.EQ.53)  WRITE  (LU,*)  FF  I  FF  9371  OR  9876 

505  WRITE  (LU  ,38'i)  DATE  ,CCOLUMN,CSAMPLE,CHOSXl,CHOSX2,CHOSX3, RES 

506  ENDIF 

507  DO  200  I»1. ,  NUMPEAK 

508  ORDER (I) -I 

509  200  CONTINUE 

510  DO  205  1=1, NUMPEAK 

511  BETA1  -  KPRIME (1,1) 

512  BETA2  -  KPRIME(2,I) 

513  BET A3  -  KPRIME(3,I) 

514  BETA 4  -  4 . 0*KPRIME ( 4 , I ) -2 . 0* (BETA1+BETA2 ) 

515  BETA 5  -  4.0*KPRIME(5,I)-2.0*  (BETA1+BETA3) 

516  BETA 6  -  4 . 0*KPRIME ( 6 , I ) -2 . 0* (BETA2+BETA3 ) 

517  BETA7  *  27 . 0*KPRIME ( 7 , 1 ) -12 . 0* ( KPRIME ( 4 , 1 ) +KPRIME ( 5 , 1 )  + 

518  1  KPRIME(6 ,1) ) +3 .0* (BETA 1+BETA 2+BETA3) 

519  X1CHOSX1/100. 

520  X2-CHOSX2/100. 

521  X3=CHOSX3/100. 

522  Y= (BETA  1 *X 1+BETA2  * X  2+BETA  3*  X3+BETA4  *X1*X  2+BETA5*X 1*X3+ 

523  1  BETA6*X  2*X3+3ETA7*X 1*X2*X3) 

514  T I  ME  ( I )  =VO  I DRT  ( 1 )  *  ( 1 +Y  ) 

525  205  CONTINUE 

526  * 

527  ******************  SELECT  ORDER  ************ 

528  * 

529  DO  215  1=1, NUMPEAK- 1 

530  DO  210  J* I +1, NUMPEAK 
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531  IP  (TIME(J) .LT.TIME(I) )  THEN 

532  ORDl-ORDER(I) 

533  ORDER (I) "ORDER (J ) 

534  ORDER(J)-ORDl 

535  Tl-TIME(I) 

536  TIME (I) “TIME (J) 

537  TIME (J) “T1 

538  ENDIP 

539  210  CONTINUE 

540  215  CONTINUE 

541  WRITE  (LU  ,386)  PL  OW 

542  DO  220  I«1,NUMPEAK 

543  J-ORDER(I) 

544  WRITE  (LU,387  )  I  ,PEAKHAME (J)  , TIMiT  { I )  .TIME  (I  )/PLOW 

545  220  CONTINUE 

546  C  CHANGE  THE  LU  NUMBERS  TO  PIT  TOUR  SYS1EM 

547  IP  (LU.EQ.6)  WRITE  (LU.381)  1  FORM  PEED  LINE  PRINTER 

546  IP  (LU.EQ.7.0R.LU.EQ.53)  WRITE  (LU,*)  PP  !  FF  9871  OR  9876 

549  380  FORMAT  (Al,  'H  ',A1, 'J ') 

550  381  FORMAT  ("1") 

551  382  PORMAT  (/,/) 

552  3e3  FORMAT  (Al) 

553  384  FORMAT  (T30,'HPLC  OPTIMIZATION' ,/,T23,15A2 ,//,T32 , 

5  54  l'DATA  ANALYSIS', /,T20, 'COLUMN  s ',A16,T45, 'SAMPLE  s  ' , A16 ,//,T25 , 

555  1 'THE  SELECTED  COMPOSITION  IS  :  ' , / ,T21 , 'MEOH '  ,T36 , 

556  l'ACN',T52,'THF',/,T21,I3,'  *',T35,I3,'  %',T51,I3,'  %',//,Tl5, 

557  1  'THE  WORST  RESOLUTION  AT  THIS  POINT  IS  ',F7.3,//) 

558  385  FORMAT  (T30,'HPLC  OPTIMIZATION ',/,T23,15A2,//,T32, 

559  l'DATA  ANALYSIS', /,T20, 'COLUMN  A16 ,T45, 'SAMPLE  :  '  .A16  ,//,T25  , 

560  \  THE  BEST  COMPOSITION  IS  : ',/ -T21 , 'MEOH ' ,T36 , 'ACH  ',T 52 , 'TUP' ,/, 

561  1T21 , 13, '  %',T35,I3,  '  »',T51,I>,'  *',//, T15,  'THE  WORST  RESOLUTION' 

562  1'  AT  THIS  POINT  IS  *,F7.3,//) 

563  386  FORMAT  (//,T27 , 'PREDICTED  ORDER  OF  ELUTION ',//,T40, 'RELATIVE  '  , 

564  1 'RETENTION  Tli-iE ',/ ,T24  ,' INGREDIENT  '  ,T40  , '1 .00  wl/mln  ', 

565  1F5.2,  '  ml/min',/) 

566  387  FORMAT  (T16 , I 2,T24 , A15 ,T42 , F6 . 2 ,T55,F6 .2 ) 

567  388  FORMAT  (//fAl,'fcdC  WORKING!  l',Al,'t@') 

568  389  FORMAT  (.U,'idC  WORKING!  ' , 1 3,A1 , *4@  ' ) 

569  RETURN 

570  END 


571  *  * 

572  **********************  SUBROUTINE  TO  PLOT  POINT  GRAPH  ****************** 

573  *  * 

574  SEMA  /TWO/, /THREE/, /FOUR/ ,/SIX/ 

575  SUBROUTINE  LCPOINT 

575  COMMON  /ONE/NUMP EAR .DATE , BESTXl, BESTX2, EESTX3 , PEAKNAME 

577  COMMON  /TWO /BEST RES .MINRES .WIDTH , KPRIME ,PEAKRT .VOIDRT 

578  COMMON  /FOUR/PEN1X ,PEN1Y , PEN2X , PEN2 Y.PEN3X ,PEN3Y ,  PEN4X , PEN4Y 

579  COMMON  /FIVE/PEN1PT , PEN 2PT,PEN3PT,PEN4PT, I SAMPLE, I COLUMN, . 

580  1  I  MEOH, I  THF , I  ACN.ISCALE 

581  COMMON  ES2,FF,L(J 

582  INTEGER  NUMPEAK ,DATE  (15)  ,  CONTRL,  I  ERR,  PEN1PT, 

583  1  PEN2PT,? EN 3PT ,PEN4PT,I SAMPLE (8) , I COLUMN ( 8 ) , I  MEOH(2),I  TH F ( 2 ) , 

584  1  I  ACN ( 2) , I SCALE ( 3) .PAGE, PLOT  "  ' 

585  REAL  WIDTH (7,35) , KPRIME (7,35) ,PEAKRT(7»35) , 

586  1  RESWAHT .MINRES (0:50,0s 50) .BESTRES, 

587  1  BEETX1, BE STX2, BEST X3,XB (0:10), YB(0:10) , XL (0:10) ,YL(0:)0) , 

588  1  XR (0 : 10 ) ,YR(0:10) , PEN1X (1 30  0) , PEN1 Y ( 1300) , PEN2X ( 1000) , 

589  1  PEN2Y ( 1000 ) ,PEN3X(1000) , PEH3Y (1000 ) , PEN4X (1000) ,PEN4Y (1000) , 

590  1  WIDE,  HEIGH'-’  ,V>1AX  ,VMIN 

591  CHARACTER  C  MEOH*4,C  ACN*4, 

592  1  C  TIIF*4  .CSAMPLE*  16  ,CCOLUMN*16  ,CSCALE*6 

593  CHARACTER  *1  ESC, FF, QUEST 

594  EQUIVALENCE  (C  MEOH, I  MEOH ) , (C  THF, I  THF),(C  ACN, I  ACN), 

595  1  (CSAMPLE, ISAMFLE)  ,  (CCOLUMN ,  I  C5LUMN)7  (CSCALETlSCALk!) 

596  *  * 

597  ****** ********************************** ******************************** 

598  *  PLOT  MINIMUM  RESOLUTION  * 

599  Gl-20. 

600  G2-350. 

601  Y_INCR“G2* (SQRT(3.0)/2)/10, 

602  DO  500  11-0,10 

603  ”«U 

504  )«INT(G1+(Y  INCR*F)/(SQRT(3.0)/2) * 

605  YB (II) -INT (Gl) 

606  XL (II) - INT (G1+ ( Y  INCR*F )/{SQRT(3.0)/l) ) 

607  Yft  ( 1 1)  -  INT  (G ).  +  (  Y— INCR*F  )  ) 

608  XR (I 1 ) -INT (Gi+ ( Y— INCR* 10+Y  INCR*F )/ (SQRT (3. 0)/l) ) 

609  YR ( I 1 ) - I NT (Gl + ( Y~ INCR* ( 10-Tl ) )) 

610  500  CONTINUE 

611  WRITE (1,590)  ESC, ESC  t  PEVERSE  VIDIO  PRINT  "PLOTfING" 

612  * 

613  ******  SELECT  PEN  COLOR  *** 

614  *  ' 

615  SCAL  E-BESTRES/ 4.0 

616  K-l 

617  L-l 

618  M-l 

619  N=1 

620  DO  510  1=0,50 

621  DO  505  J=0 , 50-1 

622  X1  =  T*  2 . /100, 

623  X2=J*  2./100 . 

624  X3  =  l  .0-X1-X2 

625  X=Gl+(1.0-X3)*G2/2+(1.0*X2*G2/2) 

626  Y=Gl+(Xl*G2*SQRr (3.0)/2) 

627  RES -MINRES ( I , J) 

628  IF  (RES  .LE. SCALE)  THEN 

629  PEN  IX ( K ) =X 

630  °P”1Y(K)*Y 

631  PEN1PT-K 
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632 

633 

634 

635 

636 

637 

638 
638 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649  505 

650  510 

651  * 

652  ***** 

653  * 

654 

655 

656  C 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680  520 

681 
682 

683 

684 

685 

686 
68  7 
688 
68  9 

690 

691 

692 


K*K+1 

ELSE  IF  (RES  .LE.  (SCALE*  2. 0) )  THEN 
PEN2X(L)-X 
PEN2 Y(L) *Y 
PEN2PT-L 
L-L+l 

ELSEIF  (RES.LE.  <SCALE*3.0))  THEN 
PEN3X(M)-X 
PEN3Y (M) "Y 
PEN3PT-M 
M-M+l 

ELS El F  (RES .LE. (SCALE*3. 75) )  THEN 
PEN4X(N)-X 
PEN4Y (N ) -Y 
PEN4PT*W 
N-N+l 

END  IF 
CONTINUE 
CONTINUE 

*  PLOT  TRIANGULAR  GRAPH  (USES  H-P  DGL  GRAPHICS  PACKAGE) 


CONTRL-O 
CALL  ZBEGN 

CHANGE  *54"  TO  FIT  YOUR  SYSTEM 

CALL  ZDINT  (54 ,CONTRL, IERR)  I  INITIALIZE  PLOTTER  ON  LU 


VM1N-0.0 

VMAX-1.0 

WIDE  »410  1  THE  MAX  X  VALUE  IN  WORLD  COORDINATES 

HEIGHT-410  I  THE  MAX  Y  VALUE  IN  WORLD  COORDINATES 

XLIt't-2  25  I  THE  MAX  X  DISTANCE  IN  MILLIMETERS 

YLIM-225  1  THE  MAX  Y  DISTANCE  IN  MILLIMETERS 

CALL  ZWIND  (VMIN,WIDE  ,VMIH , HEIGHT) 

CALL  ZDLIM  (0,XLIM,0,YLIM,IERR) 

CALL  ZASPK  (XLIM.YLIM) 

CALL  ZVIEW  (VM1N,VMAX,VM1N,VMAX) 

LMl-10/2 
DO  520  Hl-0,5 

CALL  ZMOVE  (XL(Nl)  ,YL(N1) ) 

CALL  ZD  RAW  (XB  (N 1 )  ,  YD  (tl  1 ) ) 

CALL  ZD  RAW  (XR  (111)  ,YR  (N1 ) ) 

IF  (N1.LT.LM1)  THEN 


N2-10-N1 

CALL  ZD RAW  (XL(N2) ,YL(M2 ) ) 
CALL  ZDRAW  (XB  (N  2)  ,  YB  (N2  )  ) 
CALL  ZDRAW  <XR(N2) ,YR(M2 ) ) 


ENDIF 

CALL  ZDRAW  (XL (Nl) ,YL(N1 ) ) 


CONTINUE 
CALL  ZMCUR 

CALL  ZCSIZ  (8.75,12.5) 
CALL  ZCOLR  (1) 

X=  5 


Y  =5 

CALL  ZMOVE  (X,Y) 

CALL  ZTEXT  ( 3 ,5HTHF) 
X  =  360 
7  =  5 

CALL  ZMOVE  (X,Y) 

CALL  ZTEXT  (3,3HACN) 
X=  180 


54 


Y-335 

CALL  ZMOVE  (X,Y) 

CALL  ZTEXT  (4,4HMeOfl) 

WRITE  (C  HEOH,'(I4) ')  BESTX1 
WRITE  (C*ACN, '(14)  ')  BESTX2 
WRITE  (C— THE, * ( 14) ')  BESTX3 
WRITE  (C SCALE ,  '  (F6 . 3) ')  SCALE 

XI- 30 
X2-220 
X3-240 
X4-260 
X5-280 
X6=263 
X7=300 
X8=320 
X9-310 
X10-340 

XII- 360 
X 12*  395 

CALL  ZCOLR(l) 

CALL  ZMOVE (X3,X 11) 

CALL  ZTEXT ( 18 , 18H  HPLC  OPTIMIZATION ) 

CALL  ZMOVE (X4,X 10) 

CALL  ZTEXT ( 14 , 14H  DATA  ANALYSIS) 

CALL  ZMOVE (X3,X8) 

CALL  ZCSIZ (  5. 5 , 9 . 0) 

CALL  ZTEXT (30, DATE) 

CALL  ZflOVE  (X4  ,X7) 

CALL  ZTEXT( 9 ,9H  COLUMN  :) 

CALL  ZMOVE (X9,X7) 

CALL  ZTEXT ( 15,1 COLUMN) 

CALL  ZMOVE (X4,X5) 

CALL  ZTEXT ( 9 ,9H  SAMPLE  :) 

CALL  ZMOVE (X9,X5) 

CALL  ZTEXT ( .’.5 , ISAMPLE) 

CALL  ZMOVE (X4,X4) 

CALL  ZTEXT ( 26 , 26 H  THE  BEST  COMPOSITION  IS  ) 
CALL  ZMOVE (X4,X3) 

CALL  ZTEXT (  23,2  311  MEOH  ACN  THF ) 
CALL  ZMOVE (X5,X2) 

CALL  ZTEXT (4, I  MEOH) 

CALL  ZTEXT  ( l,ltl%) 

CALL  ZMOVE (X8 , X  2) 

CALL  ZTEXT (4,  I  ACIJ) 

CALL  ZTEXT  ( 1 ,  llT% ) 

CALL  ZMOVE (XI 1 ,X  2) 

CALL  ZTEXT (4,1  THF) 

CALL  ZTEXT  (l,lff») 


********** 


PLOT  WORST  RESOLUTION  POINTS 


CALL  ZCOLR(l) 

DO  535  1=1 , PEN1PT 
WX=PEN1X ( I ) 

WY =PEN1Y ( I ) 

CALL  ZMOVE (WX,WY) 

CALL  ZMARK ( 8 ) 

CONTINUE 

CALL  ZMOVE (XI, X 11) 

CALL  ZMARK (8 ) 

CALL  ZTEXT ( 15, 15U  RESOLUTION  <  ) 
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UfJSi’XiS** -  ^  *. 

^  *-  ^  ^  v  ^  ^  •  *  -  ■  ^  ^  ~ 

754 

CALL  ZTEXT (6 , ISCALE ) 

755 

CALL  ZMCUR 

756 

CALL  ZCOLR (1) 

757 

WRITE (CSCALE, * (F6 .3) ")  SCALE*2 

758 

DO  540  I-1,PEN2PT 

759 

WX-PEN2X(I) 

760 

WY*PEN2  y ( I ) 

761 

CALL  Zf40VE  (WX  ,WY) 

762 

CALL  ZMARK  { 3) 

76  3 

540 

CONTINUE 

764 

CALL  ZMOVE (XI  ,X  10) 

765 

CALL  ZMARK (3) 

76  6 

CALL  ZTEXT(15,15H  RESOLUTION  <  ) 

767 

CALL  ZTEXT (6 , ISCALE) 

768 

CALL  ZMCUR 

769 

CALL  ZCOLR (2) 

770 

WRITE (CSCALE, * (F6 .3) *)  SCALE*3 

771 

DO  545  I=1,PEH3PT 

772 

WX*PEN3X (I ) 

773 

WY«PEN3Y( I ) 

774 

CALL  ZMOVE (WX,WY) 

775 

CALL  ZMARK (6) 

776 

545 

COMTINUE 

777 

CALL  ZMOVE (XI, X8) 

778 

CALL  ZMARK (6) 

779 

CALL  ZTEXT (15,15!!  RESOLUTION  <  ) 

780 

CALL  ZTEXT(6 .ISCALE) 

781 

CALL  ZMCUR 

782 

CALL  ZCOLR (3) 

783 

WRITE (CSCALE, '(F6.3) ')  SCALE*3.75 

784 

DO  550  1*1 ,PEH4PT 

785 

WX=PEN4X ( I ) 

786 

WY-PEN4Y(I) 

787 

CALL  ZMOVE (WX,WY) 

788 

CALL  ZMARK (9) 

789 

550 

CONTINUE 

790 

CALL  ZMOVE (X1,X7) 

791 

CALL  ZMARK (9) 

792 

CALL  ZTEXT ( 15 ,1511  RESOLUTION  <  ) 

793 

CALL  ZTEXT(6, ISCALE) 

794 

X*Gl + (1 00-DESTX3)  ‘G2/200+  (BESTX2*G2/200) 

795 

Y=Gl+ (BESTX 1*G2*SQRT (3.03/200) 

796 

CALL  ZCOLR (4) 

797 

CALL  ZMOVE (X,Y) 

798 

CALL  ZMARK (4) 

799 

CALL  ZMARK (2) 

800 

WRITE  (CSCALE, ' (F6. 3) ')  SCALE*4.0 

801 

CALL  ZMOVE (XI, X5) 

802 

CALL  ZMARK (4)  . 

803 

CALL  ZMARK (2) 

804 

CALL  ZTEXT  ( )  5 ,1511  RESOLUTION  *  ) 

805 

CALL  ZTEXT(t .ISCALE) 

806 

CALL  ZMOVE (X12,X6) 

807 

CALL  ZMARK (4) 

808 

CALL  ZMARK (2) 

809 

CALL  ZMCUR 

810 

CALL  ZD  END 

811 

CALL  ZEND 

812 

WRITE  (1,591)  ESC, ESC, ESC  !  CLEAR  THE  REVERSE  VIDIO 

813 

590 

FORMAT (A1 , *&dC  PLOTTING!  PLOTTING!  PLOTTING!  ',Al,'&@') 

814 

591 

FORMAT  (Al,'A',Al,'J#,Al,'B*) 

815 

RETURN 

*  'JCl-JSi itti "iT*. srjw Via*  --e  vu  ,-*• .  -.  v*.  --  «*•*. 


817  *  * 

818  *********************  SUBROUTINE  TO  PLOT  LINES  OF  RESOLUTION  *********** 

819  *  • 

820  $EMA  /TWO/, /THREE/, /FOUR/, /SIX/ 

821  SUBROUTINE  LCLINE 

822  COMMON  /ONE/NUMPEAK,DATE, BESTX1,BESTX2,BESTX3 .PEAKNAME 

823  COMMON  /THREE/KE SOLUTION 

824  COMMON  /FIVE/PEN1PT ,PEN2PT,PEN3PT,PEN4PT,ISAlff>LE,ICOLUMN, 

825  1  I  MEOU, I  THF , I  ACN,ISCALE 

826  COMMON  /STX/XARRAY .YARRAY ,PARRAY,QARRAY 

827  COMMON  ESC , FF ,LU 

828  INTEGER  NUMPEAK , DATE ( 1 5) ,PEN1PT,PEN2PT,P EH3PT ,PEN4PT , 

829  1  1 SAMPLE (8 ) ,1  COLUMN { 8 ) , I  MECH(2),I  THF(2),I  AGI  ( 2) ,ISCALE ( 3) , 

830  1  PAGE , PLOT, XARRAY ( 101) , YARRAY ( 101 )  7PARRAY (l5l) ,QARRAY (101) 

831  REAL  RESWANT, RESOLUTION  (7,35, 35) , 

832  1  XB (0 : 10) , YB (0 : 10 ) , XL (0:10) , YL (0 : 10 ) , XR( 0 : 10 ) ,YR(0:10) 

833  CHARACTER  PEAKNAME ( 35 ) *  15 ,CSAMPLE* 16,CCOLUMN* 16 , CSC ALE* 6 

834  CHARACTER  *1  ESC, FF, QUEST 

835  EQUIVALENCE  (C  MEOH,I  MEOI1 )  ,  (C  THF,  I  T!IF),(C  ACH,I  ACN)  , 

836  1  (CS AMPLE ,  I  SAMPLE)  ,  (CCOLUMN ,  I  COLUMN) 7 (CSCALEjl  SCALE) 

837  *  * 

838  *****************»****»*»***********.*********••*************************♦ 

839  *  PLOT  LINES  OF  RESOLUTION  * 

840  600  WRITE  (1,’C  WHAT  IS  THE  MINIMUM  RESOLUTION  DESIRED?  ") ') 

841  READ  *, RESWANT 

842  WRITE  (1,870)  RESWANT 

843  PAGE-2 

844  WRITE  (1,887)  ESC , ESC  !  CLEAR  THE  SCREEN 

845  WRITE (1,889)  DATE ,CCOLUMN , CSAMPLE, ESC, ESC  IPRINT  HEADING  CN  9876 

846  * 

847  *******  GET  BETA  VALUES 

848  * 

849  DO  750  I =1 ,NUMP EAK-1 

850  DO  745  J-I+l, NUMPEAK 

851  WRITE  (1,886)  ESC  !  TURN  OFF  GRAPHICS  DISPLAY 

852  BETA  1  -  RESOLUTION  (1,1, J) 

853  BETA2  -RESOLUTION  (2,I,J) 

854  BETA 3  -  RESOLUTION  (3,I,J) 

855  BETA4  ■  4. 0*RESOLUTION  (4,I,J)  -  2.0* (BETA 1+BETA 2 ) 

856  BETA  5  =  4. 0*RESOLUTION  (5,I,J)  -  2 . 0*  (BETA  1+BET.A  3) 

857  BETA6  -  4 . 0*RESOLUTI ON  (6,I,J)  -  2.0* (BETA2+UETA 3) 

858  BETA 7  =  2 7. 0* RESOLUTI ON ( 7 , I , J )  -  12 . 0* ( RESOLUTION ( 4 , I , J )  + 

859  1  RESOLUTION (5,1 , J) ^RESOLUTION (6,I,J))+3.0* (BETAl +BETA2+BETA3)  , 

860  * 

861  **********  CHECK  FOR  RESWANT,  -RESWANT,  AND  0 

862  * 

86  3  Y-0 

864  PLOT-O 

865  MARK  =1 

866  1 3  =  INDEX  (PEAKNAME  ( I )  ,  *  ') 

067  IF  (13. GT. 15. OR.  I3.LT.1)  13  =  15 

068  C  CLEAR  THE  SCREEN  AND  PRINT  THE  PEAK  PAIR 

869  WRITE  (1,871)  ESC  , ESC , PEAKNAME ( I )(: I  3)  .PEAKNAME (J ) 

870  WRITE  ( 1 , 8  72 ) Y 

871  GOTO  620 

872  610  IF  (PLOT. EQ . 1)  GOTO  640 

873  Y=  RESWANT 

874  WRITE  ( 1 , 872 ) Y 

875  MARK =2 

876  GOTO  620 

877  615  IF  ( PLOT . EQ. 1)  GOTO  640 
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878  Y— llBSWANT 

C79  WRITE  (1,872) Y 

880  MARK-3 

881  * 

882  *********  MAKES  PLOT  DECISION  FOR  A  PEAK  PAIR  ********** 

883  * 

884  620  DO  635  IX-0, 100,2 

885  Xl-IX/100. 

886  C1-BETA3+ (BETA1-BETA3+3ETA5-BETA5*X1) *X1 

887  C2-0ETA2-DETA  3+BETA6+(BETA4-BETA5-BETA6+BErA7- 

888  1  BETA7*Xl)*Xl 

889  C  3-0-BETA 6 -BETA7*Xl 

890  A-C3 

891,  B-C2 

892  C -Cl -Y 

893  D«B*B-4*A*C 

894  IP  (D.LT.O)  GOTO  635 

895  S-SQRT(D) 

896  Rl- (0-B+S)/ ( 2*  A) 

897  R2- (0-B-S) / ( 2* A) 

898  IF  (Rl.LT.O.OR.Rl.GT.l)  GOTO  625 

899  Q1-X1+R1 

900  IF  (Ql.GT.l)  GOTO  625 

901  PLOT-1 

902  625  IF  (R2.LT.O.OR.R2.GT.1)  GOTO  630 

903  Q2-X1+R2 

904  IF  (Q2.GT.1)  GOTO  630 

905  PLOT-1-. 

906  630  IF  (PLOT.EQ.l)  GOTO  (610,615,640)  MARK 

907  63:  CONTINUE 

908  GOTO  (610,615,640)  MARK 

909  640  IF  (PAGE .GT. 1 .AND. PLOT .EQ. l.AND.J . GT. 2 )  THEN 

910  WRITE  (1,873)  ESC, ESC, ESC, ESC  l  FORM  FEED  FOR  9876 

911  PAGE-1 

912  ELSE  IF  (PAGE.GT.3)  THEM 

913  WRITE  (1,873)  ESC, ESC, ESC, BSC  »  FORM  FEED  FOR  9876 

914  PAGE-1 

915  END IF 

916  WRITE  (53,874)  PEAKNAME ( I ) ( : 1 3) ,PEAKNAME ( J) , 

917  1  BETA  1  ,BETA2 , BETA  3 ,BETA4 ,BETA  5 ,BETA6 , BETA 7  I  OUTPUT  TO  9876 

918  WRITE  (53,891)  I  OUTPUT  TO  9876 

919  X1=0.0 

920  X2-0.0 

921  X3=l . 0 

922  Y«BGTAl*X 1+BETA2*X  2+BETA3*X  3+BETA4*Xl*X 2+BETA5*Xl*  X3+ 

923  1  BETA6  *X  2*  X3+BETA7*X1*X2*X  3 

924  WRITE  (53,875)Xl*100. , X2* 100 . ,X3* 100 . , Y  I  OUTPUT  TO  9876 

925  Xl-0.0 

926  X2-1.0 

927  X3=0.0 

928  Y»BGTAl*XlfBETA2*X2+BETA3*X3+BETA4*Xl*X2+BETA5*Xl»X3* 

929  1  BETAS  *X  2*  X3  +  BETA7*  X 1*X2*X3 

930  WRITE  ( 53  ,8 75)  XI*  100 .  , X2*  100 . ,  X3*100. , Y  I  OUTPUT  TO  9876 

931  Xl-1.0 

932  X2=0 . 0 

933  X3“0 . 0 

9  34  Y-BETAl *X 1+BETA2*X  2+BETA3*X  3+BETA4  *X 1*X2+BETA5  *X1*  X3+ 

935  1  PETA6*X2*X3+BGTA7*X1‘X2*X3 

936  WR.TE  (53,075)X1*100.,X2*100.,X3*100.,Y  I  OUTPUT  TO  9876 

937  WRITE  (53,891)  I  OUTPUT  TO  9876 

938  PAGE-PAGE+1 
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939  IF  (PLOT. HE. 1)  GOTO  740 

940  WRITE  (53,870)  »  OUTPUT  TO  9876 

941  ***  PLOT  TRIANGLE  ON  CRT 

942  Gl-10. 

943  G2-400. 

944  Y  INCR-G2* (SORT (3.0)/2>/10. 

945  DO  645  11-0,10 

946  F-Il 

947  XB(Il)-INT(Glf (Y  l!JCR*F  ) /(SQRT(  3.  0) /2 ) ) 

948  YD (II) •  JNT (G1 ) 

949  XL(Il)«INT(Glf (Y  INCR*F )/ (SQRT( 3. 0) /I ) ) 

950  YL(I1)»INT (Gif (Y"lNCR*F ) ) 

951  XR(tl)«lNT(Glf(Y"ltiCR*10+Y  I NCR*P  )/{SQRT ( 3. 0)/l) ) 

952  YR ( 1 1) -ItlT (Gif (y'INCR* ( 10-Tl) ) ) 

953  645  CONTINUE 

954  C  CLEAR  SCREEN  AND  INITIALIZE  2648  GRAPHICS 

955  WRITE  (1,877)  ESC , ESC , ESC , ESC  , 

956  LM1  *1 0/2 

957  DO  650  Nl=0,5  !  DRAW  TRIANGLE  Oti  THE  CRT 

958  WRITE  (1,870)  ESC. XL (Hi ) , YL (N 1 ) 

959  WRITE  (1,880)  ESC , XD ( Ml ) , YB (N1 ) ,ESC , XR (it 1 ) , YR(N1 ) 

960  IF  (N1.LT.LM1)  THEN 

961  N  2- 10-til 

962  WRITE  (1,880)  ESC,XL(N2),YL(N2)  ,ESC  ,XI1(N2),YD(N2), 

963  1  ESC ,XR(N2) ,YR (N2) 

964  END IF 

965  WRITE  (1,880)  ESC ,  XL  (til )  , Y L  (N 1) 

966  650  CONTINUE 

967  WRITE  (1,879)  ESC ,E SC ,ESC ,  ESC  !  LADLE  THE  TRIANGiJE 
9  68  * 

96  9  »*««•*«*  plot  LINES  OF  RESOLUTION 

970  * 

971  Y«RESWANT 

972  MARK-1 

973  GOTO  700 

974  655  Y-0 

975  MARK-2 

976  GOTO  700 

977  660  Y=- RE SWA NT 

978  MARK  =  3 

979  **“*•“*  CALCULATES  PREDICTED  RESOLUTION  VALUES 

980  700  N-l 

981  M-l 

982  DO  715  X=0,1C~ 

983  Xl-X/100 

984  Cl-DCTA  3f (DETAl-BETA 3f BETA  5 -BETA  5*Xl)*Xl 

985  C 2 -BET A 2- BETA  3f BETA6+ ( B ETA 4 -BETA 5 -BETA 6 f BET A7- 

986  1  BETA  7*X 1 ) *  XI 

987  C3-0-BETA  6 -B  ETA 7*  X 1 

988  A-C  3 

989,  B-C2 

9  90  C  -Cl  -Y 

991  D=B*B-4*A*C 

992  IF  (t.LT.0)  GOTO  715 

993  S-SQRT(D) 

994  Rl* (O-Bf S)/ ( 2* A) 

995  R2- (0-B-S)/( 2*A) 

9%  7C5  IF  (Rl  .LT.  O.Cft.Rl  .GT.l  )  GOTO  710 

997  Ql-XlfRl 

998  IF  (Ql.GT.U  GOTO  710 

999  X  2=  1-Q1 
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1000 

1001 

1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
1011 
1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 
1021 
1022 
1023 
10  24 

1025 

1026 

1027 

1028 

1029 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 
3038 

1039 

1040 

1041 

1042 

1043 

1044 
104  5 

1046 

1047 

1048 

1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 
1050 

1059 

1060 


x3-m 

XARRAY(N)-G1+ (X3+X1/2)*G2 
YARRAY (N)-G1+X1*G2*SQRT (3.0)/2 
N-N+l 

710  IP  (R2.LT.O.OR.R2.GT.1)  GOTO  715 

Q2-X1+R2 

IP  (Q2.GT.1)  GOTO  715 

X2-1-Q2 

X3-R2 

PAR RAY (M)-G1+(X3+X1/2)*G2 
QARRAY (M) -G1+(X1*G2*SQRT (3.0)/2) 

M-M+l 

715,  CONTINUE 

IP  (N.LT.3)  GOTO  72  5 
XMOVE-XARRAY(l) 

YMOVE-YARRAY(l) 

***  DRAW  RESOLUTION  LINES  ONTO  TRIANGLE 

DO  720  18-2, N-l 
XDRAW-XARRA  Y(  1 8) 

YDRAW»YARRAY(  18) 

C  DRAW  RESOLUTION  LINES 

WRITE  (1,878)  ESC,  !WOVE,YMOVE  • 

WRITE  (1,880)  ESC , XDRAW, YDRAW 
XMOVE-XDRAW 
YMOVE- YD  RAW 
720  CONTINUE 

WRITE  (1,881)  ESC ,  Y 
725  WRITE  (1,882)  ESC 

IF  (M.LT.3)  GOTO  (655,660,735)  HARK 
XMOVE-PARRAY ( 1) 

YMOVE-QARRAY (1) 

DO  730  I8»2,M-4 

XDRAW-PARRAY(IO) 

YDRAW-QARRA Y( 18) 

C  DRAW  RESOLUTION  LINES 

WRITE  (1,878)  ESC,  XMOVE,  YMOVE 
WRITE  (1,880)  ESC, XDRAW, YDRAW 
XMOVE-XDRAW 
YMOVE- YD RAW 
730  CONTINUE 

WRITE  (1,881)  ESC  »Y  I  PRINT  RESOLUTION  ON  TRIANGLE 
WRITE  (1,882)  ESC  I  TURN  OFF  GRAPHICS  TEXT 

COTO  (655,660,735)  MARK 

735  WRITE  (1,883)  ESC, ESC  !  TRANSFER  GRAPHICS  TO  9876  VIA  HPIB 

CALL  WAIT  (37,2)  l  WAIT  FOR  RASTER  DUMP  TO  FINISH 

WRITE  (1,873)  ESC, ESC, ESC, ESC  I  PORM  FEED  9876 
PAGE-1 

IF  (PLOT . EQ . 1 )  GOTO  745 

740  WRITE  (53  ,884)  PEA KNAME ( I ) ( : 1 3) ,P EARN AME (J) .RESWANT 

745  CONTINUE 
750  CONTINUE 

WRITE  (1,885) 

WRITE  (1,885)  ESC  i  TURN  OFF  GRAPHICS  DISPLAY 

READ  (1,890)  QUEST  !  DUMMY  TO  TAKE  *S"  FROM  PRINTER 

WRITE  (1,089)  ESC, ESC  I  ERASE  *S"  FROM  SCREEN 

READ  1,890)  QUEST 

IF  (QUEST. NE.  Y')  RETURN 

GOTO  600 

870  FORMAT  (//,'PEAK  PAIRS  ALWAYS  RESOLVED  GREATER  THAN  ',F5.2, 

1  '  WILL  OE  IGNORED  I  ') 

871  FORMAT  (A1  , 't!  '  ,A1  ,  M  ',//,  'FOR  PEAK  PAIR  ',Rl5,'-  ',A15) 
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1061  872  FORMAT  (//.'CHECKING  FOR  RESOLUTION  OF  *,F5.2) 

1062  873  FORMAT (M,,*Y*,A1,  'u'.M.'z'.Al,  '6p3s5dB') 

1063  874  FORMAT  (// ,T15 ,  'FOR  PEAK  PAIR  ',A15,'-  ' ,A1 5,//,T25 , '0ETA  ', 

1064  1 'COEFFICIENTS  (1-7)  ' ,/,T10 , 7F8 . 3.//.T2 3,  'MEOIl',T37,  'ACN',T52, 

1065  l'TtIF',T65,  'RESOLUTION  ') 

1066  875  FORMAT(T22,I4, '» ' ,T36 , 1  4 , , T51 , 1 4 ,  '% ',T61,F12. 3) 

1067  076  FORMAT  (//, '  ',//) 

1068  877  FORMAT  (Al ,  'll  ',A  1 ,  'J  ' ,  A1 ,  ' *m2a  0 , 0  j  lbZ  '  ,A1 ,  '*dachlnt2 ') 

1069  870  FORMAT  (A1 , '*pah ' ,214 , 'Z ' ) 

1070  879  FORMAT (\  1,  '  *d  0,0  okSTHF'  ,A1 ,  '*d  400,0  OACN  ' , A1 ,  '*(1  230,350  ' 

1071  1  'OMeOll'  ,Al,'*dlT') 

1072  880  FORMAT  (Al , ' *obh ' , 2 14 ,  '  Z  ' ) 

1C7 3  881  FORMAT  (Al , '*dsZ ' ,F 5. 2) 

1074  882  FORMAT  (Al,'*dtZ') 

1075  883  FORMAT  (Al , 'ip5D ' ,A 1 , 'fc p5u0C ' ) 

1076  884  FORMAT  ('PEAK  PAIR  ',R15,'-  ',A15,'  IS  ALWAYS  RESOLVED  ', 

1077  1 'BETTER  THAN  ',F5.2,//> 

1078  885  FORMAT  (/, 'DO  YOU  WANT  TO  PLOT  A  DIFFERENT  RESOLUTION?  ') 

1079  886  FORMAT (Al ,  '*dD ' ) 

1080  887  FORMAT  (A  1,'H',M,'J') 

1081  888  FORMAT  (T30,'HPLC  OPTIMI Z AT  ION ' ,/ ,T23, 1 5A2 ,//,T32 , 

1082  l'DATA  ANALYSIS', / ,T2 0 , 'COLUMN  :' ,A16 ,T45, 'SAMPLE  :',A16,//. 

1083  1T2,//,A1,  'll  ',A1,  '4p3s5dF') 

1084  889  FORMAT  (Al ,  A ' , Al , ' J ' ) 

1085  890  FORT-1  AT  (Al) 

1086  891  FORMAT  (T17  ,  ' - 

1087  1  ' - ') 

1088  END 
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Commander 

0.  S.  Army  Aviation  System  Command 
ATTN:  AMSAV-Q,  Mr.  E.  J.  Hollman 
St.  Louis,  MO  63120-1798 
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U.  S.  Army  Aviation  System  Command 
ATTN:  AMSAV-QE,  Mr.  J.  T.  Conroy 
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ATTN:  AMSDE-Q,  Mr.  W.  D.  Henry 
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Charabersburg,  PA  17201-4170 

Commander 
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Command 

ATTN:  AMSEL-PA,  Mr.  A.  D’Angelo 
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Commander 
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Command 
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Commander 

U.  S.  Army  Belvoir  R&D  Center 
ATTNs  STRBE-HE,  Ms.  V.  Lee 
Port  Belvolt,  VA  22060-5606 

Commander 

U.  S.  Army  Belvoir  R&D  Center 
ATTN:  STRBE-HE,  Mr.  K.  Harris 
Port  Belvoir,  VA  22060-5606 
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Dover,  NJ  07801-5000 
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Dover,  NJ  07801—5000 
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U.  S.  Army  Jefferson  Proving  Ground 
ATTN:  STEJP-TD,  Mr.  I.  K.  Peterson 
Madison,  IN  47250 
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U.  S.  Array  Yuma  Proving  Ground 
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ATTN:  STEWS-QA,  Mr.  C.  Treat 
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ATTN:  A'iSTS-Q,  Mr-  W.  G.  Creel 
4300  Goodf ellow  Boulevard 
St.  Louis,  MO  63120-5720 
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Commander 

Anniston  Army  Depot 

ATTN:  SDSAN-DQA,  Mr.  Pennington 

Anniston,  AL  36201 

Commander 

Corpus  Christl  Army  Depot 
ATTN:  SDSCC-Q,  Mr.  L.  E.  Ervin 
Corpus  Christl,  TX  78419 

Commander 

Letterkenny  Army  Depot 

ATTN:  SDSLE-Q,  Mr.  W.  J.  -.'alker 

Chambersburg,  PA  17201-4170 

Commander 

New  Cumberland  Army  Depot 

ATTN:  SDSNC-Q,  Mr.  A.  T.  Holderbach 

New  Cumberland,  PA  17070 

Commander 

Red  River  Army  Depot 

ATTN:  SDSRR-Q ,  Mr.  W.  D.  Wuertz 

Texarkana,  TX  75501 

Commander 

Lexington-Blue  Grass  Activity 
ATTN:  SDSAN-LQ ,  Mr.  P.  Jackson 
Lexington,  KY  40511 

Commander 

Sacramento  Army  Depot 

ATTN:  SDSSA-Q,  Mr.  C.  H.  Clawson 

Sacramento,  CA  95801 

Commander 

Savanna  Army  Depot 

ATTN:  SDSSV-VE,  LTC  P.E.  James 

Savanna,  IL  61074 

Commander 

Seneca  Army  Depot 

ATTN:  SDSSE-R,  Mr.  P.  W.  Chavez 

Romulus,  NY  14541 

Commander 
Sharpe  Army  Depot 

ATTN:  SDSSH-Q,  Mr.  j.  E.  Seyfried 
Lathrop,  CA  95330 
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Commander 
Sierra  Army  Depot 
ATTN:  SDSSt-QA,  Mr.  L.  B.  Lampe 
Herlong,  CA  96113 

-  -  '  ^  Commander 

Tooele  Ar.sy  Depot 

ATTN:  SDSTE-QA,  Mr.  R.  M.  Rich 

Tooele,  UT  84074 

Commander 

Pueblo  Depot  Activity 

*  '  ATTN:  SDSTE-PUG,  Mr.  P.  S.  Coen 

Pueblo,  CO  81001 

Commander 

Tobyhanna  Army  Depot 

ATTN:  SDSTO-Q,  Mr.  Vf.  J.  Lord 

Tobyhanna,  PA  18466 

Director 

U.  S.  Army  Industrial  Base 
Engineering  Activity 
ATTN:  AMXIB-MT,  Mr.  Del  Brim 
Rock  Island,  IL  61299-7260 

\  /  Commander 

Watervliet  Arsenal 

ATTN:  SMCWV-QA,  Mr.  J.  M^Iar 

Watervliet,  NY  12189-5000 

Commander 

•  Watervliet  Arsenal 

ATTN:  SMCWV-QAE-S,  Mr.  S.  Krupski 
Watervliet,  NY  12189-5000 

Commander 

Watervliet  Arsenal 

ATTN:  SMCWV-PPI ,  Mr.  P.  Phelan 

Watervliet,  NY  12189-5000 

Commander 

U.  S.  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  AMSMC-QAS-A(D) ,  Mr.  Geza  Pap 
Dover,  NJ  07801-5000 
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Commander 

U.  S.  Army  Armament,  Munitions 
and  Chemical  Command 

ATTN:  AMS MC-QAC-E ( E ) ,  Dr.  W.  J.  Maurlts 
Aberdeen  Proving  Ground,  MD  21010-5423 

Commander 

U.  3.  Army  Natick  R&D  Center 
ATTN:  STRNC-EML ,  Mr.  R.  Merullo 
Kansas  Street 
Natick,  MA  01760-5014 

Commander 

U.  S.  Army  Tank-Automotive  Command 
ATTN:  AMSTA-Q,  Mr.  L.  Barnett 
Warren,  MI  48090-2004 

Commande  r 

U.  S.  Array  Tank-Automotive  Command 
ATTN:  AMSTA-QAT,  Mr.  D.  Gamache 
Warren,  MI  48090-2004 

Commander 

U..  S.  Array  Tank-Automotive  Command 
ATTN:  AMSTA-QAT,  Mr.  F.  Braun 
Warren,  MI  48090-2004 

Commander 

U.  S.  Army  Tank-Automotive  Command 
ATTN:  AMSTA-R.CK,  Dr.  J.  Chevalier 
Warren,  MI  48090-2004 

Commander 

0.  S.  Army  Tank-Automotive  Command 
ATTN:  AMSTA-RCKM,  Mr.  D.  Ostberg 
Warren,  MI  48090-2004 

Commander 

'J.  S.  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-AD-M,  Mr.  R.  G.  Sheldon 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

U.  S.  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-AD-M,  Mr.  R.  Huddleston 

Aberdear  Proving  Ground,  MD  21005 
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Commander 

U.  S.  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-AD-M,  Mr.  W.  Deaver 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

U.  S.  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-AD-M,  Dr.  N.  Pentz 
Aberdeen  Proving  Ground,  MD  21005 

Commander 

Aberdeen  Proving  Ground 

ATTN:  STEAP-MT-G,  Mr.  C.  Klarich 

Aberdeen  Proving  Ground,  MD  21tt5 

Commander 

U.  S.  Army  Electronic  Proving  Ground 
ATTN:  STEE-PMT ,  COL  G.  L.  Kosmlder 
Fort  Huachuca,  AZ  85613 

Commander 

U.  S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-PO-Q 
Dugway,  UT  84022 

Director 

Army  Ma'.-.'ials  and  Mechanics 
Research  Center 
ATTN:  AMXMR-D 
AMXMR-DD 
AMXMR-P 
AMXMR-S 
AMXMR-0 
AMXMR-M 
AMXMR-STN 
AMXMR-ST 

AMXMR-OP,  Dr.  G.  Hagnauer 
Watertown,  MA  02171-0001 

Commander 

U.  S.,  Array  Materiel  System  Analysis 
Activity 
ATTN:  AMXSY-MP 

Aberdeen  Proving  Ground,  MD  21005 

ITT  Research  Institute 
ATTN:  GACIAC 
10  W.  35th  Street 
Chicago,  IL  60616 
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AMSMI-R,  Dr.  McCorkle  1 

Dr.  Rhoades  1 

-RK,  Dr.  Carver  3 

-RKP,  Mr.  J.  Wright  1 

-RKL,  Mr.  B.  J.  Alley  1 

-QET,  Mr.  T.  McVey  1 

-Q,  Mr.  T.  Howard  1 

-ET,  Mr.  F.  Jackson  1 

Mr.  L.  Chapman  1 

-RPR  15 

-LP  1 

-RPT  1 


Dist-9  (Dist-10  Blank) 


